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PREFACE 


This thesis entitled "Range Calculation for Heavy Ions 
in Complex Media and a Coulometric Method for Experimental 
Range Measurement" deals (a) A theoretical deduction of stopp- 
ing-pov/er equations for heavy ions in complex media for known 
molecular formula and (b) Application of coulometry in experi- 
mental measurement of ranges of fission products in copper. 

Chapter 1 explains the aims and the necessity of the 
present work in terms of the earlier related work by other 
workers . 

Chapter 2 briefly deals with the energy-loss mechanisms 
which cause an initially energetic charged particle to slow 
down inside any medium in terms of Bohr' s theory. 

Chapter 3 describes the theoretical calculation of the 
energy-loss and the ranges of heavy ions in complex media like 
nuclear emulsion. It also gives a new empirical method of cal- 
culating the mean ionization potential for any elemental medixim. 

In Chapter 4 the various general experimental procedures 
employed for the range determination of heavy ions and fission 
products are described. 

Chapter 5 presents the coulometric method used in the 
present work for the differential range measurement of the 
fission products ®^Sr, ^^Zr and in metallic copper. 
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CHAPTER 1 


INTRODUCTION 


Theoretical calculation of the energy- loss of energetic ; 

charged particles in matter has been one of the important unsolved | 

problems in the realm of physics for a long period of time. In 

1 

1913, Niels Bohr first formulated the method of solving it by a 

2 

purely classical method and was followed by Bethe who used a 

3 

quantum mechanical approach. In 1948, Bohr looked afresh at 

the problem in detail and gave a solution which was consistent 

1 2 

with his earlier treatment as well as with Bethe 's quantum 
mechanical formulation. However, the stopping-power equation of 

3 

Bohr contained some summations, the evaluation of which could 
not be done by any rigorous theoretical procedure. In an approxi- 
mate manner these are obtainable from a knowledge of the velocity- 
distribution among the orbital electrons in an atom. In the case 

3 

of a heavy ion (Z> 40) penetrating a heavy medium (2> 40), Bohr 
himself suggested the use of the Thomas-Fermi statistical approac'*' 



2 


for this purpose. For comparatively lighter ions penetrating 
light media, no statistical approach is suitable so that even 
a simple approximation is theoretically ruled out. Mukherji 

4 

and Srivastava, using a semi-empirical approach, showed a way 

3 

of making Bohr' s equation practically useful in the case of 
partially- stripped heavy ions like fission products and similar 

5 

accelerated heavy ions. Later they tried to extend it to the 

case of heavy ions which are almost completely stripped because 

of their very high velocity. Although, in general, the stopping- 

5 

power equations deduced by these authors predicted ranges and 
energy- loss of heavy ions in elemental media, which were in good 
agreement with the corresponding experimental values, there were 
some notable exceptions. The present work, in part, deals with 
an investigation of the reasons of the disagreement between the 
calculated and the e^qjerimental values in the case of some ele- 
mental media. After locating the reasons and making suitable 
modification in the computational procedure, the stopping-power 
equations have been used for the calculation of ranges of heavy 
ions in complex media like mylar, polyethylene, tissue material 
and nuclear emulsion. Although nuclear emulsions have been used 
for charged-particle identification from the lengths of the 
tracks, there has been no appropriate, even semi- empirical, gene- 
ral equations relating the track lengths to the velocity of a 
charged particle of known mass and atomic numbers. The available 
track-length versus energy equation of Heckman ^ is based 

on many dubious and arbitrary assumptions. The most serious 



3 


drawback of the equation is that it is of an "ad hoc" character. 

The method used in this work has been step-wise. The equations 

4 5 

developed by Mukherji and Srivastava and Srivastava and Mukherji 
for stopping-power calculations were first modified so that there 
was good general agreement between the experimental energy- loss 
in elemental media measured by other authors and those calculated 
by these equations/ without any exceptio n. In the next step, 
energy- loss in relatively sirrple compound media like mylar and 
polyethylene were calculated, using Bragg's additivity rule. This 
rule assumes that the heavy ion sees no chemical binding in the 
compound medium and the stopping-power of the compound medium is 
the sum of the stopping-powers of the constituent elements, 
weigtfed according to the nximber of atoms of each consituent per 
molecule. The ranges of heavy ions in these media have been 
computed from the stopping-powers obtained in the aforesaid way. I 

I 

After good agreement was obtained between the experimental and | 

i 

calculated values in those cases, ranges in nuclear emulsion 
were computed using two different approaches, both leading to 
practically the same final values. ; 

The other part of the present work consists of testing a 
new experimental technique for measuring the ranges or depth of 
penetration of fission products or of similar heavy ions in copper. 
Of the two general methods of range determination (i.e., "integralj 
and "differential") the differential method requires removal of 
very thin layers of the "catcher" material and subsequent deter- 
mination of the concentration of the heavy ion in the layers by 



4 


suitable methods like counting of radioactivity. Although it is 
well-known that metals like copper and silver dissolve quantita- 
tively when used as an anode in an electrolytic cell, no attempts 
have been made to apply this principle for differential range 
measurement. In this work, the conditions for dissolving uni- 
formly thin layers of copper has been investigated. The agree- 
ment between the experimental and calculated ranges has been 
quite good, indicating that the electrolytic method may possibly 
be applied in the case of other metals similar to copper. 

The results of the present theoretical and experimental 
investigations have already been reported in a journal.* 


* 1. Nucl. Instm. and Methods, 159 , 421 (1979) . 

2. Nucl. Instm. and Methods (accepted) (as part of the 

Proceedings of the 8th International Conference on Atomic 
Collisions in Solids) . 



CHAPTER 2 

THEORETICAL ASPECTS OF STOPPING PROCESS 

2.1 ENERGY LOSS PROCESSES 

The mechanism of energy loss of an energetic charged 

particle involves four different types of collisions with matter 

8 

through which it passes. 

2.1.1 Inelastic collision with atomic electrons 

For swiftly moving charged particles, inelastic collisions 
with bound atomic electrons of the medium being penetrated is 
the most important energy loss process. The atomic electrons 
experience a transition to an excited state (excitation) or to 
an unbound state (ionization) . 

2.1.2 Elastic collision with the atomic nucleus 

In this type of interaction, also called nuclear collision, 
the incident particle gets scattered in the screened field of the 




6 


nucleus of the stopping atom. During the process the kinetic 
energy and the momentum are conserved. Nuclear collisions become 

Q 

predominant at low ion- velocities, 2, 2x10 cm/sec. 

2.1.3 Inelastic collision with the nucleus 

In a close non-capture encounter with the nucleus of the 
stopping atom the incident charged particle experiences a deflec- 
tion. In some, but not all such deflections, a quantum of radia- 
tion is emitted and a corresponding amount of kinetic energy is 
lost by the incident charged particle. Since the cross-section 
for this type of emission is inversely proportional to the mass 
of the charged particle, this process is not important in the 

g 

case of heavy charged particles. 

2.1.4 Elastic collision with the atomic electrons 

The incident charged particle can suffer an elastic colli- 
sion with the electrons of the struck atom. Both kinetic energy 
and the momentum are conserved, and the energy transfer is less 
than the lowest ionization potential of the bound electrons so 
that the interaction is really with the atom as a whole. This 
type collision is important only in the case of very low energy 
electrons as projectiles. 

All the four energy- loss processes described above can 
contribute to the energy loss suffered by the incident particle. 
Their relative contributions, however, depend upon the energy and 
the type ( heavy or light) of the charged particle. For heavy 
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charged particles, inelastic collisions v/ith the atomic electrons 
of the stopping medium, also called electronic collisions, domi- 
nate at high ion velocities, while at low velocities, of the 
order of the most loosely bound atomic electrons of the stopping 
medium, only nuclear collisions are important. 


2.2 ELECTRONIC STOPPING PR0CESS5S 

A brief account of the various electronic stopping equa- 
tions is given in this section. 


2.2.1 Bohr's classical stopping power formula 

A detailed discussion of the interaction between a moving 

charged particle and an electron of the stopping medium was given 
1 

by Bohr in 1913. Using nonrelativistic classical mechanics and 
Rutherford's model of the atom, Bohr obtained the following 
equation for ^ , the energy loss per unit path length due to 
electronic collisions; 


dX 


4Tr 

mv^ 


n 



1.123 mV" 


ze 


0 ). 


. .. ( 2 . 1 ) 


where ze is the ionic charge of the incident particle, is the 
atomic number of the medium, n is the nximber of atoms per unit 
volume of the stopping medium, m and e are the mass and charge 
of the electron respectively, V is the ion velocity, and lo is 
the cyclic frequency of the s orbital electron of the stopping 
atcxn. The sum includes all the orbital electrons of the stopping 
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atom. Eq. (2,1) can also be written as 


dE 4Trz^e^ _ „ ^ 1.123 mV' 

_ = ^ n Z 2 In — ^ T 


( 2 . 2 ) 


mV 


ze (13 


where o) is the geometric mean cyclic frequency of the orbital 
electrons of the stopping atom, given by 


Z 2 In (j) 



V 

4^ 

s=l 


In (0 

s 


(2.3) 


2.2.2 Bethe's stopping power formula 

The first purely quantvim mechanical treatment of energy 

loss of charged particles due to electronic collisions was given 
2 

by Bethe based on the Born approximation. The use of the Born 
approximation recguires that 


ze 


•hv 


« 1 


,(2.4) 


where ze and V are the charge and the velocity of the charged 

particle. This condition is usually satisfied at high velocities 

and small charges of the moving particle. Under these conditions 

2 

the non-relativistic stopping -power formula is given by 


dE _ 4Tr z^e^ „ V ir, 2 mV^ 

dX „2 n ^2 I 

mV 


(2.5) 


where I is the mean excitation potential for the atoms of the 
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stopping medium. For relativistic velocities of the charged 
particle the energy-loss is given by 


dX 


4 


TT 

mV^ 



{In 


2 mV^ 



where P = / and c is the velocity of light. 


( 2 . 6 ) 


2.2.3 Bloch ' s equation 

The condition for the validity of the Born approximation, 
^2 2 

« 1 , and the classical approximation, comple- 

3 

mentary to each other. A bridge between Bohr's equation and 

9 

Bethe's equation was provided by Bloch who took into considera- 
tion the perturbation of wave functions the atomic electrons of 
the stopping atom in the presence of a charged particle. The non- 

9 

relativistic stopping -power formula due to Bloch is 


g = iJLzV ^ ^2 {In Rtf, (1+i ^)} (2.7) 

mV 


where is the logarithmic derivative of the gamma function and 

2 - 
Rt|, denotes the real part of if, . For » 1 Sq. (2.7) reduces' 


to Eq. (2,2) and for 


ze 

fiV 


« 1 it reduces to Bethe ' s formula. 


Eq. (2.5). For intermediate values of , Eq. (2.7) is approxi- 

mately identical with Eq. (2.2) or Eq. (2,5). 


Bloch^ made another interesting contribution by showing 


that for atoms which can be adequately described using the 
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Thomas-Fermi statistical model, the mean excitation potential I 
is proportional to the atomic number 

I = k (2.8) 

The "Bloch constant" k has to be determined experimentally. 

2.2.4 Bohr's semi-quantTom mechanical treatment of electronic 
collisions 

The classical equation of Bohr and the quantum mechanical 
equation of Bethe are based on the conditions that x >> 1 and 
X « 1 / respectively, were x is given by 

V 

X = 2 z — (2.9) 

In Eq. (2.9) z is the charge of the ion in units of the electronic 

2 

charge e, V is the ion velocity, and ^ . That the two equa- 

tions predict the same energy-loss for x = 1 is, in a sense, acci- 
dental, since they are derived for values of X which are mutually 
exclusive^ 

3 

In 1948 Bohr attempted a more detailed analysis of the 
problem using a semi-classical approach. A brief derivation of 
these stopping power equations is given below: 

When a swift charged particle is penetrating matter the 
loss in kinetic energy will be primarily due to collisions with 
the atomic electrons of the penetrated medium. The collisions 
can be divided into two types, namely, "free collisions" and 
"resonance collisions". 
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Free collisions - When the time of interaction between the 
charged particle and the electron of the stopping atom is very 
small compared to the period of oscillation of the electron in 
the atom, to a first approximation, the influence of the atomic 
binding forces can be neglected and the electron can be considered 
as free. Such collisions are called "free collisions;' Free colli- 
sions are accompanied by ionization and excitation of the mediTm. 

Resonance collisions - When the time of interaction between the 
charged particle and an electron is very large compared to the 
period of oscillation of the electron in the atom, the field 
exerted by the moving charged particle will be practically uniform 
over the entire atomic region. This type of collisions are called 
"resonance collisions". 

"til. 

The ionization potential of the s electron, I , is 

s 

given by 


I 


s 



( 2 . 10 ) 


where U is the velocity of the s^^‘ electron of the stopping 

s 

medium, and m is the electron mass. If the length a^ is defined 
by 


a 


s 


h 

mU 

s 


( 2 . 11 ) 


then a is a measure of the radius of the orbit of the s 
s 

electron. In quantum mechanical terminology, mU^ can be 
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3 

considered as the expectation value of the moraentum of the 
electron and as the accuracy with which the electiron can be 
localized without the uncertainty in its energy exceeding I . 

If 0 ) is the cyclic frequency of the s^ electron, the ioniza- 

O 

tion potential, I , can also be written as 

3 

w, (2.12) 

As pointed out earlier, the problem of electronic collision 
effects produced by fast moving particles depends essentially on 
the value of the quantity X . Hence, it is convenient to treat 
the two cases x > 1 and X <1 separately. 


X > 1 

The classical orbital picture of the electrons can be used 

I 

to describe the individual encounters between the particle and 

3 

the atomic electron. The collision time, t , for coulombic 
interaction for such systems is given by x =p/V, vrhere p is the 
impact parameter. For this to be of the order of the period of 
oscillation of the atomic electron, the impact parameter p=d 
is given by 


P = 



(2.13) 


where m is cyclic frequency of the s electron. 

3 

The energy transfer, T, to a free electron by a charged 
particle for an impact parameter p is obtained from the relation 
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T 


^24 
2 z e 


mV 


1 



(2.14) 


If p = i is the impact parameter for which T = I , the ioniza- 
s s 

tion potential of the electron. 


i 


s 



(2,15) 


where the collision diameter, b, which is the shortest possible 
distance of approach in a head-on collision between the electron 
and the incident particle is given by 

b = ^-^4“ . (2.16) 

mV 

When V >> U , d >> a and a > b, d can be taken as the effective 
^ s s 

adiabatic limit, and for impact parameters greater than d the 
energy- loss will be negligibly small, and as a first approxima- 
tion can be neglected. 

Using Eqs. (2.10), (2.11) and (2.12) d^ can be v/ritten as 


d 


s 





where a new parameter n defined by 

& 



(2.17) 


(2.18) 
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has been introduced. Using Eqs. (2.15) and (2,17) one obtains 





-1 

X 


(2.19) 


For ri„> Xr io always less than d and the collisions with 
impact parameters p <ig may be considered as free and the eneirgy 
transfer T will be always greater than the ionization potential, 
I . The energy-loss due to free collisions is given by 


T 

(AE ) =n.AX.B. ln{(^)} 
s f e I 


( 2 . 20 ) 


In Eq, (2.20) is the maximxim energy transfer to the electron 

2 

and is equal to 2 mV , n is the number of atoms per unit voltame 
of the stopping medium. Axis the distance traversed by the 
charged particle, and Bg is defined by 


Be = 


2 ir z^e'^ 
mV^ 


( 2 . 21 ) 


Substituting the value of and ^ Eq. 
as 

2 

(Ae^)^ = n*AX*B_ In } 

s f ^s 

2 

= n* AX - B In n 
e s 


(2.20) can be written 


( 2 , 22 ) 


For i > a , the encounters with P> i„ are of a simple resonance 
s s s 

type. The contribution to energy loss due to resonance collisions 
is given by 



15 


I 

(A Eg) = n.AX .Bg In , (2.23) 

s 

where D is the energy transfer corresponding to impact parameter 

dg. Substituting the values of and from Eq. (2.10) and 

(2.14), and expressing the logarithiTiic arguments in terms of n 

s 

and X one obtains 

(A Eg)^ = n.A X. Bg In ( X~^) • (2.24) 

If X in which case i > d , the values of p for which 

S S3 

the collisions can be considered ‘free’ no longer extend to i^ 

since for p > d the duration of the encounter exceeds 1/ u .On 
^ s s 

the otherhand, all collisions for which p> d will not be of a 

s 

3 

purely adiabatic type. According to Bohr the limiting value 

d* of the impact parameter for which the probability of ioniza- 
s 

tion can still be assimed to be unity is given by 

(d * ) ^ = i .d . (2,25) 

s s s 

If D” is the energy transfer corresponding to the impact para- 
s 

meter p / the energy-loss due to free collisions is given by 

(AE^) =n.AX.B_ In (T / D" ) 

= n. AX. Bg In ( x“^) , 


(2.26) 
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holding for T^» D^* or x « The total energy- loss to the 

til 

s orbital electron E due to free collisions and resonance 

o 

effects will be the sum of Eqs. (2.22), (2.24) and (2,26). 




AE 


n»Ax.Bci In 




In Ti 


rx 

fe] ] 


(2.27) 


where the quantity inside the square bracket, if less than unity, 
should be replaced by unity. 


X <1 

For X < 1 the orbital picture of classical mechanics 
completely fails in accounting for individual collision effects. 
However, for X <1 the separation between the "free" and 

3 

"resonance" collisions is possible to a first approximation, 
and these correspond to p < a and p > a respectively. The energy- 
loss due to free collisions is still given by Eq. (2.20) while 
that for resonance effects is given by 

S 

A 

( Ae ) = n • Ax • B_ In , (2.28) 

Sr ^s 

where A' is the energy transfer calculated according to classical 
s 

mechanics for p = a . On substituting the values of A' and D , 

s s s 

and simplification, one obtains 


( AE ) = n . A X • In 

Sr ^ 


2 

^ * 


(2.29) 
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Thus/ for x<l^ the total energy-loss, E , to the s^^ orbital 

3 

electron of the stopping medium is given by 


A = ( AE ) + (AE ) 

s S f S 

= 2 n . A X. B_ In n ^ . 

^ s 


(2.30) 


A comprehensive 
Eqs. (2.27) and 


3 

energy-loss equation can be written by combining 
(2.30) : 


Ae 

Ax 


Il« B 


1 

S lnTig^[3Q"^+ 2 In {tg^r^T } . (2.31) 

s s s 


where the terms within the square brackets if less than unity 
are to be replaced by unity, and the summation is to be carried 
over all orbital electrons which interact with the moving charged 
particle. (^) / the rate of energy-loss per unit path length, is 
called the stopping power of the raedi'um for the charged particle 
concerned. 

2.2.5 LSS THEORY 

10 

In 1963 Lindhard et worked out a theory for energy- 

loss of swiftly moving charged particles based on the Thomas-Fermi 
model of the atom and the theory of quasielastic collisions between 
the charged particle and the atoms of the stopping medium. Accord- 
ing to this theoiry the electronic stopping-cross section, per 

0 

^ 10 . 

atom IS "given by 
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Q 2 

8 ire a 


^2 




-v_ 

V. 


(2.32) 


bolding for the velocities of the charged particle, v < V^. 
In Eg. (2.32), and are the atomic numbers of the charged 
particle and the stopping medium respectively, a^ and are the 
radius and the velocity of the electron in the first Bohr orbit 
of the hydrogen atom, e is the electronic charge, and 5 is of 
the order of 


2,3 NUCLEAR COLLISIONS 

The coulomb interaction between the screened nuclear 
charges of the incident charged particle and the atoms of the 
stopping medium, often referred to as "nuclear collision" 
becomes the dominant mechanism of energy-loss only when the 
velocity of the charged particle becomes less than the velocity 
of the most loosely bound electrons of the stopping .atom. It 
is simply Rutherford scattering between two charged nuclei modi- 
fied by proper screening by charge clouds of the orbital elec- 
trons . 


In terms of energy transfer, T, the Rutherford differen- 
tial cross-section da for collision between two particles carry- 
ing charges q^e and q_2^ is given by 


da (E,T) 


2 2 4 

2 TT qi q2 e 


m2 V 




wq. ^ q.2 e ^ 

^2 ' 


E 


(2.33) 
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where E is the kinetic energy of the incident particle in the 
laboratory co-ordinates/ and are the masses of the inci- 

dent and struck particles respectively. In the case of a colli- 
sion between a low-velocity heavy ion and an atom of the stopping 
medium the presence of the bound atomic electrons causes shield- 
ing of the nuclear charge and hence the interaction potential 
has to be suitably modified. The exponential screening poten- 

3 

tial due to Bohr between two atoms with nuclear charges Z^e 
and Z 2 e is given by 


V(r) 



(2.34) 


where 'a ‘ is the screening parameter defined by 

3 


^s = 


O 




(2.35) 


The applicability of Rutherford's equation depends on the ratio 

r between *a ' and the collision diameter *b' for an unscreened 
s 

field: C = b/a. If C « 1 and x » 1 / the interaction v/ill be 

3 

confined mostly to the unscreened part of the field. For 
strong screening, i.e., ? » 1, the collisions may be considered 

as occurring between two rigid spheres, the radii of the spheres 
will, however, depend upon the energy of the incident particle.- 
Taking these factors into consideration and using Rutherford' s 
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3 

scattering theory# Bohr derived an expression for the rate of 
energy-loss due to nuclear collisions: 


dX 


2 2-4 

2 irZ^ Z 2 ® 


n 


ln( 


M1M2 

Ml +M2 




), 


( 2 . 36 ) 


where a is the screening parameter defined earlier, 
s 



CHAPTER 3 


■ CALCULATION OP HEAVY ION RANGES IN 
NUCLEAR EMULSION AND OTHER CCMPLEX MEDIA 

3.1 INTRODUCTION 

The stopping-power equations, one due to Bohr^ and the other 

2 • 3 11 

due to Bethe, have their restricted regions of validity. 

Recently, starting with the classical stopping-power forrmila due 
3 5 

to Bohr, Srivastava and Mukherji deduced a set of stopping- 
power equations for calculating the electronic stopping-powers 
of heavy ions at all ion velocities greater than e . 'While 
these equations predicted stopping-powers and ranges quite well 
for fission fragments and heavy ions in light stopping media 
(i.e., atomic numbers less than 14) , considerable disagreements 
were observed in the case of heavier media. This necessitated 
a re investigation of these stopping-power equations. 
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3.2 DERIVATION OF STOPPING POWER EQUATIONS 

The energy-loss per -unit path length dE/dX for a heavy 
ion of ionic charge ze and velocity v due to electronic colli— 

3 

s ion has been given by Bohr : 


dE 

dX 


2 4 

2 TT z^e n 
inV^ 


{ £ In 
s 


{%^ ra) 




(3.1) 


where n is the nxirriber' of atoms of the medium per tin it volume, 
e and m are the charge and the mass of the electron, x and n_. 
are defined by 


X = 2 z 


Vo 

V 


(3.2) 


_ 2V 


(3.3) 


where U„ is the velocity of the s 
s 


th 


orbital electron of the 


stopping atom, and V^=e'^/h. In Eqn. (3.1) the quantities within 
the square brackets, if less than 1, should be replaced by unity.' 
Using this condition and the definitions of x ^nd , .the second 
ternT in Eq. (3.1) can be split into two terms as follovvB: 


S ince. 


2L 


u_ 


2V X 


-1 


for U < 2V x”^/ -^ < 1 and for > 2V 
s Hg s 




U =2V X- 


, o r y ® 

S '' S 1_T1 I 


In Ti, 


U3=o 


s s 


U^=2V X 
s 


> 1. Hence, 


-1 


ln(ii. 


X"^) (3.4) 
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Therefore Eq. (3.1) can be written as 

_2 U =2V x“^ 

ln(n ^CxH )+ 2 In Y)l + 

U3=0 s 

u =u " 

^ ® 3-11 

I In (r,; X ^)} 

U^=2VX ® 

S 

(3.5) 

I tl 3 

where U and u are the upper cut-off velocities '' vfiich render 

o o 

the logarithmic terms zero conforming to the fact that physically 

no negative energy-loss is possible. For convenience Eq. (3,5) 
will be written as 


dE 

dX 


2 tt z\-‘^n 
mV^ 


{ 


Uo=U., 
s s 


U3=0 


^ - lJ ., g (J + J + J ) 

mV i z J 


dX 


where 


s s 


-2 


= 2 InCug^UxI 1 / 


U =0 
s 


Jr-. = 


U3=2V X 


U = 0 

S 


in n 3 " , 


® ® 3—1 

2 In (n^ X , 

\J^=2V X ® 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


and, further, the sum of the three terms and will be 

represented by S, 
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S =* + J 3 . (3.10) 

O 

In Eq. (3.5) the summation terms, as Bohr'^ showed, could be 
replaced by integrals if the velocity distribution fxmction for 
the orbital electrons were known. As a first approximation, for 
atoms which can be adequately described using the Thomas -Fermi 
statistical model, Bohr gave the following formula for n(Ug) : 


n(U^) 



(3.11) 


where n(U ) is the nuntoer of orbital electrons with velocity less 

o 

2/3 

than u for V<U<ZV . Zis the atomic number of the 
s o s ^ o 

atom. For light atoms, vhich are not expected to obey the Thomas- 
Fermi statistical model Eq. (3.11) is not valid. Various empiri- 
cal estimates for n(u ) have been made.^ In the present work an 
empirical formula 


n(Ug) = f(Z) (3.12) 

with f(Z) = 0,3634 has been xised (cf. Section 3.5). 

In order to obtain a set of general stopping-power equa- 
tions it is necessary to evaluate the sxmmation terms correspond- 
ing to different physical conditions determined by X , z, v and 
the atomic number of the stopping medium, Z 2 . The two cases 
X >1 and X <1 are treated separately for convenience. 
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X > 1 


i) ^ I ^2^0 X 

can be written as 

J. = S In in j x~h = £ In ( 2V/U x ) ^ (3,13) 

u^=0 ® u=o ® 

S 


since the logarithmic term remains positive for all values of 

-1 

U from U =0 to U„=2V X = Z^V , the summation between the 

S o o ^ O 

limits U = 0 to U.=u’ is equivalent to the inclusion of all the 

o S S 

Z 2 electrons of the stopping atom. The velocity of the K-shell 

electrons of the medium of atomic n\arriber Z^ is taken as Z^V^. 

z z o 

Since n„ can be written as 


2V 

U 


2 mV 


_ 2 mV 


(3.14) 


where I is the ionization potential of the orbital electron, 

s 

for one obtains 

s=Z2 2 

J-. = £ In ^ 

0=1 T X^ 


= Z 2 in 


2 mV 
I X' 


(3.15) 
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where I is the geometric mean ionization potential of the atoms 
of the stopping medium, and is defined by the equation 

s=Z2 

Z, In I' = z In I . (3.16) 

s=l ^ 

Since V '^- 2^0 ^ been set as the upper limit, for one 

obtains 


J 


2 


Ug=2V X 

S In n ^ 

U =0 ® 

s 


£ In (2V/U ) 

U =0 ® 

s 


I In (2mV^/I ) 
s=l 


Z 2 In (2m vVl ) . 


(3.17) 


Since the lower limit of in corresponds to the maximum 
possible value of U , becomes redundant. Therefore 

o O 


S = Ji + Jo + J-, = 2 


^2 


(2mV^ / i X) 


(3.18) 


ii) 


ZoV X 
2 O 


> V> 


(Z„-2) V^X 

Z O 

2£(Z2) 


In this case the ion may not be able to ionize all the 
atomic electrons of the meditim. In evaluating the contribution 
of term it is useful to replace the suunmation by integration 
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Using the differential form of n(u ), viz,, 

s 


dU 


ancu^) - fCZj) — 


(3.19) 


Since the two K-shell electrons of the atom are practically im- 


. .12 


shielded Mukherji suggested that their contribution to the 


stopping power can be computed separately. Accordingly, can 


be written as 


(Z2-2)V^ 


f(Z2) 


1/3 2 1/3 3 

J, = ; In (2V/LT X ) an(U ) + 21n ( 2V/Z,V X ) 

O s s ^ o 

2V X { 3 . 20) 


The upper limit of the first integral stands for the velocity 
th 


of the ( 22 ~ 2 ) orbital electron of the medium as given by 
Eq, (3.12) and the second teinn represents the contribution of 


.5 


the two K-shell electrons. Earlier, Srivastava and Mixkherji” 

^2^0^ 

evaluated in the ion velocity region — ^ > V > ( 2 , ) — 


without considering the possibility of the occurrence of the 


-1 <V2>''o 

condition 2V X ^ 


f(Z2) 


. (Z.,-2)V^ 

If 2V X > 


£(Z2> 


the, inteara] 


in Eq. ( 3 . 20 ) becomes zero, and further, 
its 
^2^0 ^ 


if V < 


z o 




itself becomes zero. Thus the following two cases arise if 


(a) V> 


> V > - 

22'^oX 


(Z„-2) 
z o 


2 f(Z2) * 


In this case 


j- = 2 In { 2V/Z.,V^ ^ 

o Z O 


( 3 . 21 ) 
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As far as is concerned since the upper cut-off value uj equals 
-1 -1 

2V X / and as 2V X is greater than the velocity of the 
th tHo-t of 

(Z2~2) electron but less^han the K-shell electrons, the simma- 
tion would include only the outer {^^-2) electrons. Thiis, 

u’=2v x“^ 

J-1 = Z on (n/ X 


s — ^ 

= s In (2invV x^) 
s=l ® 

= (Z 2 - 2 ) In (2rtiV^/ X^) (3,22) 


where is the geometric mean ionization potential of the outer 
(Z 2 “ 2 ) electrons of the stopping medixim and is given by the rela- 
tion 

s=Z2“2 

(Z^-2) In i, = Z In I . ' (3.23) 

^ ^ s=l s 


The value of has been obtained from Kukherjis formula 


12 


Z,-2 

I, = 13.6 { ^ > 

2.717 f(Z2) 


(3.24) 


In a similar manner one obtains for J 2 

,-l 


U =2V X 
S 


s=Z2-2 


= S In = Z In ( 2mV / 1_ ) 

2 u =0 ^ s=l ® 

s 


= ( Z 2 - 2 ) In ( 2mvV \ ) 


(3.25) 
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Therefore, 


,2 


S = J^+ J^+ J 3 = 2 {Z^~ 2 ) In (2mvV x) + 


2 In ( 2 V/Z 2 V^ ^ 


(3.26) 


ZoV^ X‘ 
2 o 


.1/3 


(b) V ^ 2 • Iri this case becomes zero, -while and 


J 2 are given by Eqs. (3.22) and (3,25) respectively. Thus, 


2 , - 


S = J 3 _+ J 2 = 2 (Z 2 - 2 ) In (2inV / x) 


(3.27) 


(iii) 2 f(Z 2 ) 


2 o 


1/3 


tin 

since 2V X is less than the velocity of the (Z 2 - 2 ) 


electron, bo-th and J 2 -would require an upper cut-off value 


of 2V X . Therefore, for , J 2 and one gets 


U =2V X 
s 


-1 


J-, = S 

11=0 

s 


-1 


U =2V X 
s 


-1 


In (2V X VUg) =2 J- 




ln(2v X Vu ) dn(U ] 

-'0 C i 


U =0 
s 


2f(Zj V^VX 


-1 




4f(Z^) X ^ V 


V. 


J In (2V X /u ) dU = 

s s 


V 


u^=o 


(3.28) 


U =2V X 


-1 


U =2V X 
s 


Jn = 


£ In (2V/U ) =2 J 


In { 2 V/U ) dn (U ) 

o o 


,-l 


2f(2,) “s=2V X 

f In ( ^ ) dUg= 


4 f(Z 2 ) X (1+ In X) 


V 


U3=0 


V. 


V 


(3.29) 
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U_ = 


( Z,-2) 

I o 


s ~fTzp~ 


^3 = 


2V X 


-1 ® 


“•o 

Z O 


Sfiz^) V 

_ 

o 


^ ^2-^> 


2V X 


; in (-^^) au, + 21n ( 

-1 U xl/3 X s ' xl/3 / 


3f(zJ (Z„-2)V^ 

I r I O 


-1 


v_ 


f(Zo) 


- 2V X + 2V X In X 


-1 ^-2/3 


( Z„-2) V ( Z^-2) V X 

f(Z2) 


1/3 


2V f(Z2) 


+ 6 In i~~) - 2 In X> 

^2^0 


(3.30) 


.?hus 


2f(Z ) V 

S = J2+ J2+ J3 = 3(Z2-2) +3(Z2-2) In ( z -2) V 

2 O 


2V ^2^ ^ 

+ 61n 2 V V X ~ ^2 

2 o o 


(3.31) 


Z^V X^/^ (Z„-2)V^ X^/^ 

(i^) _2_o_^ > V > 2 ° 


2f(Z2) 


and J2 would be still given by Eqs. (3.28) and (3.,2^ 
respectively. As far as is concerned, since V < 


Z2V0 x^/-' 
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the second term on the right hand side in Eq. (3,20) becomes 


zero. Therefore, for one obtains 

f ( ZJ 

3 

Jo = I In ) dn(u ) 

—1 TT ® 

2V X ^ 


,iV!lZo 


2V X ^ + 2V In X^/2 


(Z,-2)V„ (Z,-2)V„ 

^ q , 2 o , 

f(Z 2 ) 2 v f(Z 2 ) ^ 


(3.32) 


Hence, 


J]^+ J 2 + J 3 = 3(Z2~2) + 


2f(Z2) V 


(Z2-2)V X 

-3(Z2-2) in t 


C3.33) 


1/3 

<3^) V < ' X 

In this case and J 2 are given by Eqs, (3.28) and (3.29) 
respectively, requires an upper cut-off value of uj'= 2V ^ 

Thus, for one obtains 


3f(Z2) 


2V X 


-1/3 


■^3 V 


° 2V X~^ 


^ "^ 73 ^ = 


-4f(Z„)V X' 


^s 


-In X 


3f (ZJ 


— (2V -2V X“^) 


(3.34) 
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and hence 


2 = '^ 1 + >^ 2 + '^3 = 


V 


(3 X 


-1/3^ 


x”^) 


(3.35) 


X < 1 

For X < 1 

u„=u ' 

s s 

= £ In and (3.36) 

U ^=0 

u =2V x“^ 

J 2 = 2 In (3 '3) 

u =0 ® 

S 

since Ug = 2v X ^ is greater than the maximum allowable value of 
Ug' = 2v X , for ^ < 1/ becomes redundant. 


In this case even the F-shell electrons of the medium can 
be ionized. Therefore, the suiimation in and Jae over 

all the Z 2 atomic electrons. 

2 2 V ^ ^~^2 2 

Jl = J 2 = £ In (^) =2 £ ^ In (2inVVls) 

S=1 *^3 S=1 

= 2Z2ln (2mvVi) (3.37) 


Therefore, 

S = J 2 = 4 Z 2 In (2mvVi) (3.38) 

where I has been defined alreac^ by Eq. (3.16) . 
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( i — £—9 . > V > ^ , P 

2 V > 2f(Z2) 

__ 

Since V < the maxirnmn velocity imparted by the 

particle is less than the orbital velocity of the F-shell elec- 
trons. Therefore,- the F-shell electrons can not participate in 

the energy loss process. However, the velocity of the (Z^— 2)^^ 
(Z2-2 )v^ ■ 2 

electron is „ \ — . Therefore all the electrons upto the 


th 


f(Z2) 


( Z2“2) electron would be ionized. Hence, 

2 


s=Z2-2 

jj = = S in (|2) 

s=l s 


s~Z2“2 

= 2 E In 
s=l 


2mv' 


= 2 (Z 2 - 2 ) In 


2mV^ 


(3.39) 


where has been already defined by Eg. (3.23) 


2mv 


S = J^+ J 2 = 4(Z2-2) In ~ 


(3.40) 


(iii) V < 


(Zj-av^ 

zfCZj) 


This condition requires an upper cut-off value of 2V for U 


U =2V 
s 


u^ = o 


Jl = J 2 - ^ in (f^) 

s 


2f(Zj 


2V 


X in (Si) an(Uj 


V. 


4f (Z 2 ) V 
V 


■u 


(3.41) 


S = J 2 = 


8f(Z2) V 
_ 


(3.42) 
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3.3 HYDROGEN AS THE STOPPING MEDIUM 

The general stopping power equations derived above for 
multielectron atoms have to be modified for hydrogen as the 
stopping rnediim. Since hydrogen atom has only one orbital elec 
tron the summation in Sq. (3.1) becomes redundant and U has to 
be replaced by V^. Thus, 

S = ln [f j' } 

1_ g-J. 

= in Djh (3.43) 

'^o o — 

where/ as before/ the (Quantity within the square bracket, if 
less than 1, should be replaced by unity. The two cases X > 1 
and X < 1 will be treated separately. 


X > 1 


In this case S is given by 


S=ln j . 

c c * — 


(3.44) 


Depending upon the magnitude of the ion velocity two cases arise; 


(i) V> 




V Y 
o ^ 


under this condition — g — < 1 and therefore. 


, ,2V^ , - o ir, 2mV‘ 

S = In (~) + In (“) - 2 In -- 

''o o 


(3,45) 
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where I is the ionization potential of the hydrogen atom. 

V 


V X 

( ii) — > V> ^ 


V Y 

o^ 


For — ^ — >V the first term in Eq. (3.43) will be zero. 

V X V X 

Further/ for V'<— ^ > 1 and/ hence/ the square brackets 


in the second term can be removed. Therefore/ for S one gets 


S = 3 In y., . (3.46) 

X 

o 

238 

Since even for heavy ions like U-ion x = 8 the lower limit 
V 

Y = ... 9 — . .^ (3.46) will hold for range calculations 

for any ion. Below V nuclear collisions accompanied by 

large angle scattering dominate the energy-loss process and hence/ 
as far as range calculation is concerned/ the contribution of 
electronic collision towards energy-loss below V = can be 
neglected. 


x< 1 

In this case both the logarithmic terras become identical 
and S is given by 


S = 2 in ( 

o 

= 2 In ( ) (3.47) ' 
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3.4 stopping power EOTTATTOktc; 

After substitution of the expressions for S from Eqs. 
(3.18), (3.26), (3.27), (3.31) , (3.33) , (3.35), (3.38), (3.40), 
(3.42), (3.44), (3.45), (3.46) and (3.47) in Eq. (3,1) and 
conversion to proper units and simplification one obtains the 
following expressions for the stopping-powers, in Ifev-cm^/mg, 
valid under the conditions specified. 


(i) X>1, V>,, 


Vo 


^ = 63 . 65 11^39 

dX . „2 ^2 ^°%0 


A 2 V 


I X 


(3.48) 


ZV (Z„-2)V^X 2.V Z^V 

(ii)X>l, -1^ X>v>-2_^ : 


-2 2 


^ = .63.65 z j- 11.39 V“ . ^ , 2V 

c3X ^^^2 ^^^2 2) log^o + 3 log^^ ^ 


ZoV X 

z o 


1/3 


(3.49)1 


(2o-2)V X Z^V_ X^/^ 


(iii) x>l, -3.^ >^Vo 


2f(Z2) 


2-0 “ 

2 > ^ ^ “ 2f(Z^) 


dE _ 63.65 - 11.39 




( Zj-2) log;^o ^ 


X 


(3,50); 


Z^V X 


1/3 


x>i, :z:2 


, 'V2)VoX ZgV^X (V2)V„X 

^ o-p/''7 ^ ' 3 > V > 


2f{2,^) 


2 ' 2f(Z2) ■ 2 

Eq. (3.5 0) will also be applicable in this region. 
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V-., ,^2^0 , <V2>V,X (22-2)V„X __ xV" 

' 2 ^ 2f ( z^ — “Tfr^r” — 2 — 


i = -r :| ^ »<22-2) + 3(22-2) 1" T2-32TV 


2f{Z^ V 


A^v 


-f- 6 In 


f ( Z^ V 

2 -y + 2 y X ^2 ^ ^ 

2 O O 


2V 


<3.5X) 


(vi) 


Z 2 V (Z.-2)V^X Z \f xl/3 (Z,-2)V^ X^^^ 

X>l,-2_o_ <__2_^,_2o 2 


2f(Z2) 


2f{Z2) 


dE 

dX 


= ^3.79 z' 


2f(Zj V 

{3(Z2-2) + 

/-V ^ 


3 (Z 2 - 2 ) In 


(Z„-2) (V X' 
2V f ( Z 2 ) 


,1/3 


(3.52) 


( vii) X > 1, 


Vo (Z,-2)V„X (V^>V . (“2-2) V 


1/3 


> 


•'2 o 
2f(Z2) 


2f (Z^) > 2fTz^ 


Eq. (3.52) would be also applicable in this region. 


(Z„-2)V^ X^/2 

(viii) X> 1, V < 2f("Z2) ' 

i = ^(^2> (3 .X-3/3 . x-3, 

(ix) X < 1, V > ; 

2 

dE _ 63.65 2 Z 2 ^ 11.39 

ax j 


(3.53) 


(3.54) 
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(x) 


X <1, 


Z^V 
2 o 


> V > 


(Z,-2)V^ 
I O 

2f(Z2) 


dE 

dX 


63.65 


(Zj-a 


l°9l0 


11.39 

T 

"■1 


( 5 . 55 ) 


(xi) 


X <1, 


V 


< 


(Z.-2) 

2 O 

2f(Z2) 


dE 

dX 


50 , 6 

A 2 V 


2 ^ f ( Z 2 ) 


( 5 . 56 ) 


The following equations hold in the case of hydrogen as 
the stopping medium; 


(xii) X > 1, v> 




dE ^ 63.65 z^ 11.39 


(3.57) 


^o ^ ^o 

(xiii) X>1, -%-> V> -2-2 


1/3 


dE _ 47.74 11.39 

dX ~ ::2 2/3 


V 


I X 


(3.58) 


V. 


(xiv) X < 1/ V>” : 


dE 

dX 


63.65 z^ 11.39 V^. 

-2 log^O I 


(3.59) 
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In all the above equations V and (= 2.185 x 10® cm/sec) 

are expressed in mits of 10® cm/sec, and I, and I (= 13.6 eV) 
in eV. A 2 and are the mass and atomic nimbers of the stopp- 
ing medium, z is the ionic charge in units of the electronic 
charge e. The proper equations for I, , f(Z 2 )and z are given 
in next sections (3.5) and (3.6). 


3.5 MEAN IONIZATION POTENTIAL 

In Bethe's stopping power equation, already mentioned 
earlier (section 2.3.2), 


dE 

dX 


4-ir z^e^ 
m 


Z In 


2 mV 


(3.60) 


I stands for the mean excitation potential of the medium of 
atomic nvimber Z. This parameter is defined by 


Z In I 


= 2 
i 


f. In (h 


0 ) • ) 
X 


(3.61) 


where the summation is extended over the various virtual osci- 
llators of strength f^^ and frequency These oscillators 

are attributed to the transition probabilities of the individual 

atomic electrons c It is quite difficialt to obtain accurate 

11 

estimates of I on purely theoretical grounds. Most workers 
have to obtain the values of I for a given medium from Bethe’s 
equation with the help of experimental values of the energy-loss 
of a known ion with known initial velocity. Bloch^ put forward 
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the idea that, on the basis of the Thomas-Permi statistical 
model, the mean excitation potential should be proportional to 
the atomic number 

I = KZ (3.62) 


where K is a constant. The difficulty, however, was that when 

the values of I for many media were obtained from Bethe' s 

(3.P 

equation using experimental values of , the values of K were 

found to vary with the media. Further, the variation of K with 

8 

Z is such that it is difficult to make any interpolations or 
extrapolations on the basis of known values of K to obtain the 
mean excitation potential for an imknown medium. 


1 2 

Mukherji obtained a relatively simple formula for the 
calculation of I. He used the assimmption,v initially made by 

•a 

Bohr, that, as a first appreximation, the mean excitation 
potential I is identical with the mean ionization potential I, 


s=Z 

Z In I = Z In I = Z In I , (3.63) 

s=l ® 

th 

where I is the mean ionization potential of the s electron 
s 

and the summation includes all the Z electrons of the atom of 
the medium. The s\ammation was replaced by an integration 
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Z In I = 2 In ° (In 1 ) dn(U ) 

s=l ® U =0 s s 

s 

U =ZV 
s o 

= / In 4 m ) dn (U ) , (3.64) 

U =0 ^ s s 

s 

where m is the mass of the electron^ U is the velocity of the 

3 

th 

s orbital electron and n(U ) is the niomber of electrons having 

s 

velocity less than a given velocity U . For n(U ) Mukherji and 

s s 

4 

Srivastava have given the following expression mentioned earlier. 


n(U^) = f(Z) ^ 

where f(Z) = 0.28 Z^^^ for Z 4 45.5 . 

f(Z) = Z^'^^ for Z >45.5 . 


(3.12) 

(3.65) 

(3.66) 


However, Eqs. (3.65) and (3.6C) do not hold for either the inner- 
most or the outermost electrons of the atom. After appropriate 

- 12 

modification, the expression for I was given as 

2 

Z In I = (Z-2) In [l3.6 { TTfiTTTTP ^ ^ 

(3.67) 

where I is given in eV. The values of I predicted by Eg. 3.67) 
were in fair agreement with the values of I obtained throu^ 
experiments in the cases of light and heavy elements. There 
was considerable disagreement in the case of medixam- weight 
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elements of atomic number, 13<Z<74. Significantly, the calcu- 
lated ranges of heavy ions in nickel and silver as stopping 
media show large deviation from the experimental values. Since 
all stopping-power equations contain f(Z) directly or throu^ I 
and it was assumed that a new prescription is needed for the 
calculation of f(Z) . There are three reasons for the assumption. 
Firstly, there is a large difference between the experimental 

values of the mean ionization potentials and those calculated 

12 

using Mukherji’s equation in the case of elements of atomic 
number Z in the region 74>Z>26. Secondly, from the theoretical 
point of view Eqs. (3.65) and (3,66) were obtained from a conside- 
ration of the interaction between a moving ion and only a few 
outer orbital electrons, except in the case of very light media, 
and, hence, should be valid only for the outermost electrons, 
whereas their use in Eq. (3,67) implies their validity for all 
orbital electrons except the two K-shell electrons. Lastly, 
since f(z) largely represents the "effective quantum ntmiber" of 
the electrons, which, as Bohr argued wouild be unity for the 
outermost and the innermost orbital electrons, the values given 
by Eqs. (3,65) and (3.66) for f(Z) are possibly overestimations 
for medium-weight and heavy-weight elaiients, particularly at 
high ion energies when a majority of orbital electrons are 
involved in the interaction. Thus, one should take values of 
f(Z) averaged over all orbital electrons; these woudd be smaller 
than the values obtained tising Eqs. (3.65) and (3.66), would lead 
to larger values of I, and would, therefore, be consistent with 
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1 O 

the experimental values. To obtain the new expression for 
f ( z) Eg. (3.67) was assumed to be basically valid and that by 
substituting the experimental values of the mean excitation 
potentials in place of I in that equation one woiHd obtain* the 
averaged or new values of f(Z) . The experirttental values taken 
from Bakker and Segre,^^ as calculated by Bethe and Ashkin^^ and 
from the compilation of Ba’rkas and Berger^^ are shown in Table 
3.1. Figure 3,1 shows a logarithmic plot of the new f(Z) values 
against log Z. A least squares fit yields the relation 


f(Z) = 0.3634 


(3.68) 


which on siibstitution in Eg. (3.67) / gives 

2 

Z In I = (Z-2) In { ^ ^ 2 in (13.6 Z^) 

(3.69) 

and 


For Z<13 and Z>74, the values of f(Z) from Eg, (3.68) differ 

very little from those obtained by the earlier equations Eqs. 

(3.65) and (3.66), but for 13 <Z<74, there is a B-17% lowering 

in the values of f(Z). Table 3.2 lists the newly calculated 

values of I for some elements along with the corresponding 

1 6 

accurate experimental values from Andersen ^ and Siz5rensen 

17 

and TUidersen which have not been used as inputs in obtaining 
the relation given by Eg. (3.68). Since the earlier expressions 
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Table 3.1 

The experimental mean excitation potentials ( in eV) xised 
in obtaining the expression f(Z)= 0.3634 


Element 

Atomic 
nvurriber, Z 



Li 

3 

34 

38 

Be 

4 

60.4 

60 

c 

6 

76.4 

78 

N 

7 


85 

0 

8 


89 

Me 

10 


131 

Al 

13 


163 

Ar 

18 


210 

Pe 

26 

241 

285 

Cu 

29 

27 6 

314 

Ag 

47 

418 

487 

Sn 

50 

463 

516 

W 

74 

655 

748 

Au 

79 


797 

Pb 

82 

705 

826 

U 

92 

811 



^Bakker and segre (ref. 14). 
Barkas and Berger (ref. 15). 



log f(Z) 


-0.4 


- 0.6 


Fig. 3.1 Plot of log,o 

the data fra 
Berger .(Ref 
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Table 3,2 

Comparison between the calculated values of I with the 
experimental values of mean excitation potential I 


Element 

Atomic niimber, Z 

I (ev) 

Experimental 

I (ev) 

Calculated 

Ca 

20 

195.0® 

231.2 

Sc 

21 

217.0® 

241.0 

Ti 

22 

229.0® 

250.7 

V 

23 

238.7® 

260.2 

Cr 

24 

259.1^ 

269.8 

Mn 

25 

273.8® 

279,2 

Fe 

26 

280.8® 

288.6 

Co 

27 

299.3® 

297.9 

Ni 

28 

303.4® 

307,1 

Cu 

29 

321.7® 

316.3 

Zn 

30 

323.4^ 

325 . 4 

Ag 

47 

466.8^ 

474.0 

Au 

79 

750.0^ 

733.8 

Pb 

82 

766.0^ 

757.4 

U 

92 

831.0^ 

835.0 


. Andersen, H. Sorensen and 

P. vajda (ref. 

16) . 


Sorensen and H.H. Andersen 

(ref. 17) . 
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for f(Z) were deduced from the ranges of fission products, the 
new f(Z) values from Eq. (3.68) would yield considerably larger 
values for calculated fission product ranges. As fission 
products are essentially heavy ions having low initial kinetic 
energies in the range 0.5 -l.O iXfeV/amu, in all stopping power 
calculations Eq, (3.68) has been used at ion energies greater 
than 1.0 MeV/amu while below that energy Eqs. (3.65) and (3.66) 
have been used. 

3.6 CALCULATION OF THE ION CHARGE 

At low ion velocities, conparable to those of outer atomic 
electrons, capture and loss of electrons by the penetrating ion 
complicates the energy-loss calculations. With increasing atomic 
number of the ion, over wide ion-velocity region, charge fluctua- 
tions become important. However, the processes involved in 
capture and loss of electrons by a moving ion are in general so 
intricate as to defy any precise and comprehensive formulation. 

According to Bohr, “ an ion penetrating a material medium 
loses all electrons whose orbital velocities, U , are less than 
the ion velocity V. Therefore, the charge ntomber of the ion ( z) 
will be the same as n(U ) with U V* in Eq. (3.12) 

® o o 

z = n(Ug=v) = ^ (3.71) 

o 

where \ is the atomic number of the ion and is given by 

Eq, (3.68) with Z=Zj^, Equation (3.71) has been found to hold 
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good only for ions which are less than half stripped, that is, 

4, 18 

z < , For higher ion velocities at vhich the ion is more 

than half stripped, Srivastava and I'^ukherjf have given the 
following expression for calculating the ion charge, 

z = Z^{l-c exp (-2f( Z^) V/Z^V^) } . (3.72) 


Since the ion charge can be predicted quite well usina Eq.(3.71) 

2 


ZiV 

for ion velocities V 4 xT-f ir / that is vhen z < ~ , the cons- 


2f(Z) 

tant c in Eg. (3,7 2) can be obtained by assuming that at 
Z^V 

V = 2fTz) ' (3.71) makes a smooth transition to Sq. (3.72) 


= {1 - 2.039 exp (-2V f ( Z^) /Z^V^) } . 


(3.73) 


Though semi-errpirical in nature, the good agreement between the 
calculated and the experimental values justifies the use of 
Eqs. (3.71) and (3.73) from the point of view of practical 
utility. 


3.7 RANGE COMPUTATIONS 

The energy transferred to the atomic electrons of the 
stopping medium in a single collision by a fast moving heavy 
ion is quite large compared to chemical bond energies. Therefore, 
the fast moving ion does not essentially see any chemical bond 
between the atoms of the coir^xDund, Hence, the stopping power of 
the compound medium can be calculated by summing up the stoppiftg 
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powers of the component elements taking into accoimt the partial 
density of each elemental component as follows. The stopping 


power 


dE 
dX/i. 






of the i atomic species in its pure elemental 


state at a particular ion~energy s is civen bv 



_ 2 IT 2^6^ 


mV 


n. 


(3.74) 


where n . is the number of atoms per unit volume and s • is defined 

JL 

by Eg. (3.l0) , Eg. (3,74) can also be written as 


fdE\ 1 


2Trz^e'^ 

mV^ 


P.N 

— 1_ e 
Aj_ ^i " 


or 


1 1 

f dE] 

7 

_ 2Trz^e^ 

N 

i <»=/ 


mv^ 

^-i 

where ^ 

■i is 

the 

density of 

the ; 


. th 


atomic mass , and N is the Avogadro number. For the complex 


medium the stopping-power 


L(i)i 


■-lE 


(3.75) 

atomic species of 
ber. For the comi 
at ion energy E is given by 



2£z^ef 

mV^ 


n 


Y. 3 . 
1 1 


(3.76) 


where n is the number of molecules of the compoimd per unit 
c 

th 

volxme/ and Yi is the number of atoms of tbe i ‘ atomic species 
per molecule. If A^ is the molecular mass and is the density 

of the complex medium. Eg. (3.76) can be rewritten as 
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dE \ 

dxjc 


JE 


2Tr2^e^ 

mv^ 


PcN 


- £ Y S . 
c i ^ 



A~ 


£ 

i 


A.Yi 


[(■ 


n 


i5B\ 
dX ■ 


;aJ.E 


(3.77) 


after substituting for from Sq. (3,74) . The range R of the 
ion can be computed numerically; 


5 



the energy interval 6E should be sufficiently small (‘^O.Ol MeV) 
so that the stopping power remains sensibly constant over this 
small interval; E^ is the initial energy of the ion and E^ is 
the energy corresponding to ion -velocity . After the ion has 

slowed down to the velocity its energy loss occurs largely 
through nuclear collisions. These involve random large-angle 
scatterings, the net result of v^ich is that the ion does not 
penetrate much in it's original direction after the velocity 
becomes less than . The numerical computations were done 
using an IBM 7044 computer. The computer was provided with the 
equations (Eqs. (3.4^) - (3.59) ) for stopping powers, Eqs. (3.71) 
and (3.73) for calculating the ionic charges, along with the 
velocity regions in v^iich they are valid. I and were calcu- 
lated using Eqs. (3.69) and (3.70 respectively. The f(Z) values 
in all the above equations were obtained by mean of Eq. (3,68) . 
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The energy of the ion was changed in steps of 6S=0,0l i-teV from 
the initial energy (= lo Mev/amu) to the lower limit 3^ 
corresponding to the velocity . The stopping powers were 
obtained using appropriate equations for each energy interval 
and the ranges were calculated using Eq. (3.78). Ranges were 
also calculated with 6E =0.001 i-teV in a few cases. The ranges 
so obtained, however, differed from those calculated with 6S = 
0.01 MeV by < 1%. 


3.8 RESULTS AND DISCUSSION 


3 ..8 . 1 Ranges in Elemental ^dla 

Since many of the stopping power equations presented here 

5 

differ from those obtained by Srivastava and Mukherji, it was 
necessary to check the validity of these equations for simple 
elemental media before using them in the case of complex media. 
The energy-losses for ^^0, and ions with 

initial energies of 10 MeV/amu in passing through various thick- 
ness of oxygen gas were calculated using the procedure described 

in the earlier section. The calculated energy-loss curves aire 

19 

compared with the experimental values from Roll and Steigert 
in Pigs. 3.2 and 3.3. In Fig. 3.4 the calculated energy-loss 
curves for '^*^Ar and ions in nitrogen and argon gas are 


compared with the corresponding experimental values from Martin 


20 


and Northcliffe.^"^ Nitrogen, oxygen and argon are typical light 
media and, as the Fias. 3.2, 3.3 and 3.4 shc^J 40BM^ent 

CfcHTHAl I 

A 6 6 81S 


No. 
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200 r 



Thickness of Oxygen gas (mg/cm^) 

Fig. 3.2 Energy-loss curves for and ions 

in Oxygen gas. Experimental data ore from Roll 

and Steigcrt (Ref. 19) 
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Thickness of Oxygen gas {mg/crn^) 

Fig. 3.3 Energy-loss curves for’^0, and ions in 
Oxygen gas. The experimental data are Iron Roil 

and Steigert (Ref. 19) 
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between the calcxilated and the experimental values is quite good. 

For light stopping media/ however/ even the earlier calculations^ 

were in good agreement with the experimental values. It vras in 

the case of heavy stopping media like nickel that the earlier 

results were in considerable disagreement with the experimental 

values. Hence/ it is crucial to ascertain the applicability of 

the present stopping power equations to heavy stopping media. In 

Pigs. 3,5, 3.6 and 3.7 the calculated energy loss curves for 
11 „ 12 ^ 14 ., 16 ^ 19 ,, , 20 

B, C/ N, 0/ F and Ne-ions in nickel are compared with 

. IQ 

the experimental data from Roll and Steigert. The calculated 

^oss curves for C- ions in silver using both the old and new 

f(Z) values are shown in Figure 3,8 along with the experimental 

21 

values from Walton and Hubbard. The good agreement between the 
calculated and experimental values for nickel and silver shows 
that the present stopping power equations can be used for medium- 
weight elements as stopping media. 

3.8.2 Compound Media 

22 

Schambra rt al. have measured the energy-loss suffered 
by accelerated heavy ions of initial energies of lo Jiev/amu in 
passing through different thicknesses of mylar (<^0^8 *^4^ 
polyethylene, (CH 2 )j^. The energies of the incident and the emer- 
gent iqn were measured with magnetic spectrograph. The experi- 

12 16 

mental energy-loss curves for C, O- ions in polyethylene/ 
and ^^0, ^°Ne- ions in mylar are compared with the calculated 

ones in Figs. 3.9 and 3.10 respectively. For calculating the 





Ion energy (McV) 



3.6 Energy-loss curves for and ions in Nicke 
The experimental data are from Roll and Steigert 
(Ref. 19) 




^'^nerav (MeV) 
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Thickness: of Nickel (mg/cm^) 



Energy- loss curves for and ions in Nickel. 

The experimenta! date ore from Roll and Steiger! (Ref. 19) 



Ion energy {Mev /a mu) 
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Fig. 3.8 Range-energy curve for ions in silver. The experimental 
data are from Walton and Hubbard (Ref. 21). The solid and 
dashed lines represents the ranges calculated using the new 
and old f(z) values respectively. 



I 
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Polyethylene thickness (mg/cm^) 

Fig. 3-9 Energy-loss curves for’^C and ^^0 ions In '^c'yet^^ylene. 

The experimental data are from Schambra etal (Ref. 22) 



Ion energy {MeV/amu) 


20n \16o 


Mylar thickness (mg/cm^) 

Fig. 3.10 Energy-ioss curves for ’^0 and ^®Ne ions in Mylar 
The experimental data are frexn Schambra etal (Ref. 22) 




62 


energy—loss/ the coinpovind media was treated as a homogeneous 

cif the constituent elements and chemical binding was 
neglected. Energy—loss curves have also been calculated for 
simulated tissue material having the chemical composition 

^7^70^^2*^32 ^ ions and ^*^Ne and compared 

with the corresponding experimental values from Brustad^^ in 
Fig. 3.11. In all the cases the agreement between the calculated 
and the experimental values appears to be good. ^Vlar, and poly- 
ethylene and the simulated tissue material, however, contain only 
light elements and the present results are not significantly 
different from those obtained by srivastava and I'^herji^ for 
these particular media. 

3,8.3 Nuclear Research Emulsions 

Nuclear research emulsions contain finely divided (micron 

f 25 

sized) ^ins of silver bromide uniformly dispersed in gelatin. 

Since the energy-loss suffered by an ion in passing through a 
single silver bromide grain can be quite large, the nuclear emul- 
sions possibly can not be regarded strictly as homogeneous media 
at molecular level. Ranges in G-5 nuclear emulsion were calcu- 
lated by using two different procedures. In the first procedure 
the nuclear emulsion was considered as a homogeneous medium viiile 
in the second procedure the heavy ion was considered as passing 
alternately through segments of silver bromide grains and gelatin 
each retaining its own bulk density. 



Stopping Flower CMc^V an^/mrp 



Fig. 3.11 Pio^s of stopping power of '“8, ‘®0 end ions 

against their penetration depth in tissue material. The 
experimental data arc from Brustad (Ref. 24) 
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Prom the composition of the "standard" G-5 nuclear emul- 

. . 

s ion, given in percentage by weight, , one obtains the following 

/ 

empirical foimnula (Appendix) for G-5 nuclear onulsion of density 

3 

3.845 g/cm ; 


^0.994 ^0.006 


'1.370 ^"0.314 ’'^0.013 ^3.185 °0.932 


when the emulsion is under vacuum during bombardment ^ due to the 
evaporation of water, the density changes to 4.0 g/cm^. Let 

3 

xg/cm be the weight of water lost vhen kept under vacuum. Since 

each gram of water removed causes a reduction in volume of 

3 25 

0,73 cm , as obtained from the dehumidification curves, one 

obtains 


3.815 - X 
1 - 0.73X 


4.0 


(3,79) 


and hence x = 0,0962. Taking into account the evaporation of 
water, the new empirical formula is 

^'^^^0.995 ^0.006 ^,370 ^0,314 '^0,013 ^2,545 %,62S 

Assiiming the empirical formula of gelatin to be C 3 _^ 37 q ^o. 314 
^O 013 ^2 545 °0.626 ' neglecting the small atomic propor- 

tion of iodine, the path lengths through gelatin and silver 

25 

bromide were obtained following the procedure of Barkas. 
average grain diameter <D> of silver bromide in G-5 nuclear 


The 
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25 

emulsion is 0.275Mm, with a standard deviation o= 0.0022Mm, 

Assuming random traversal of grains the average path length <6 > 

”* s 

through silver bromide is given by^^ 

<6g> = I <D> = 0.1833 ;«m, (3.80) 


If n be the number of grains the ion passes through when it 
traverses unit distance, the average distance t between the 
centres of two successive silver bromide grains is given by 

t = Vn . (3.81) 


The value of n has been calculated from the relation 


25 


n 





<D>‘ 


) (1 - 


3 O' 


<dX 


(3.82 


where C is the fraction of the voliime occupied by silver bromide. 

In calculating C the change in volume due to evaporation of water 

during bombardment under vacu'um has to be taken into account. 

The weight of water lost per cm'' of the original emulsion is 

0,0962g, and the corresponding decrease in volume obtained from 

3 3 

dehumidification curve is 0,0702 cm . Therefore 1 cm of the 

3 

oricfinal emulsion v/ill occupy a volume of 0.9298 cm during 
bombardment. Thus, if the change in density of the emulsion 
under vacuiam is attributed to thee^faporation of water only, 

0.9298 cm^ of the emulsion will contain 3.1407 g of silver bromide. 
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Assuming that silver bromide retains its bulk density of S.AIZ 
3 

g/cm auring bombardment, for C .one obtains a vlue of 0,5239. 

Hence, L is equal to 0.3712 Mm and the average path length 
through gelatin, t is given by 

<6g> = L- > = (0.3721 - 0.1833) M m 

= 0.1888m m. (3.83) 

I ! 
s 

I 

The heavy ion passes alternately throu^ an average distance of 
0.1888 Mm of gelatin and 0.1833 Mm of silver bromide. ' 

For computational purpose eacli silver bromide and gelatin 
segment was theoretically subdivided into several sections so 

I 

that the energy-loss suffered by the ion in passing through eadi j 
section did not exceed O.Ol MeV. The range was obtained by noting J 
the number of segments of silver bromide and gelatin through v^idh i 

the ion passes before its velocity is reduced to V . The ranges S 

O I 

calculated using this procedure for the ions | 

^^0, and ^*^Ne are shown in Pigs. 3,12, 3.13 and 3,14.along- 

26 

with the corresponding experimental values from Roll and steigeut. 
These calculated ranges did not differ sensibly from those cal- ; 

f 

culated by the first procedure stated earlier. The calculated ; 

ranges are systematically larger than the experimental values ! 

by vA3Mm at ion energies below t/'4.5 ^5eV/OTU for all ions other ' 

than ^®Ne for which there is almost a constant difference of i 

5 M m at all energies. The small disagreCTiait at low energies ; 

i 

can be due to : (i) the distortion of gelatin during the ! 



Range (/Jfn) 
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I 



Fig. 3.12 Range -energy curves for and '^0 ions in G-5 
emulsion. The experimental date are from Roll 

Steigert (Ref. 26) 


nucleor 

and 



Range (jjm) 
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Ion energy (M 'V/amu ) 


Fig. 3.13 Range-energy curves for and ions in G-5 

nuclear emulsion. The experimental data are from 
Roil and Steiger! (Ref. 26) 
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development process involving htimidification and heating viiic±i 
can introduce some error in the measured track lengths, parti- 
cularly for low energy ions, (ii) for heavy ions the uncertainty 
in their charge state, and, hence, in their kinetic energy, 
increases as the ion energy decreases below <^4 MeV/amu. 

The overall agreement between the calculated and the expe- 
rimental energy-loss curves appears to be good for all kinds of 
stopping media, and hence, one can conclude that these stopping 
power equations along with the modified formula for ionization 
potential are of general applicability. 



CHAPTER 4 


TECHNIQUES OP RANGE MEASUREMENTS 


4.1 RANGE CONCEPTS 

A swiftly moving charged particle loses kinetic energy by 
various types of interactions which have been discussed earlier# 
and comes to a halt after penetrating a certain distance of the 
stopping medium. The penetration depth is called the range. 
However# since the stopping processes are statistical in nature# 
there is no unique range for a large number of particles having 
the same initial kinetic energy and the same initial direction 
of penetration. Hence# it is customary to define ranges in 
terms of the ''average" range# the "most probable" range and the 
"median" range. The most probable range is given by the peak in 
the range distribution curve, and the median range corresponds 
to the penetration depth beyond which 50% of the incident parti- 
cles penetrate. If the range distribution curve is gaussian# 
then all these ranges will have the same magnitude. 
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It should be noted that experiments give the projection 
of the actual path in the initial direction of penetration. Only 
for those cases where slowing down of the energetic particle 
tabes place via a large nxrriber of small-angle scattering events 
does the range and the actual path length have almost the same 
magnitude. For heavy charged particles having high kinetic 
energies/ the energy-loss is predominantly due to inelastic 
collisions with the electrons of the stopping medium and the 
above approximation holds good. However/ low energy charged 
particles in penetrating heavy stopping media lose kinetic energy 
by collision with atoms of the stopping medium which is accompa- 
nied by large-angle scattering and/ hence, for these ions there 
is a considerable difference between the actual path length and 
the range; since theoretical calculations refer to the total 
path length/ the experimentally obtained ranges may require 
corrections before any valid comparisons can be made between the 

two. 


4.2 RECOIL TECHNIQUES 

When the target nuclei, e.g,, U, undergo fission, the 
fraction of a particular fission product escaping from the target 
depends upon the direction of emission, its range, as well as the 
thickness of the target. Depending upon whether the target thick- 
ness is negligibly small or very large compared to the range of 
the fission product of interest (also called recoil product) , the 
different experimental methods for range measurement can be 
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roughly divided into two categories, namely, thin target methods 
and thich target methods, respectively,"^ 

4.2,1 Thin Target Methods 

In this method the target has a thickness much less than 
the range of the recoil product of interest and hence the emerg- 
ing recoils can be assumed to preserve the original range dis- 
tribution and the angular distribution. 

(i) Thin target-thin catcher differential method- A well- 
collimated beam of recoil products emerging from a very thin 
target is allowed to penetrate the catcher material at a known 
angle of incidence. The catcher can be a stack of thin foils or 
it can be a thick foil from which layers of uniform thickness 
can be removed. When a stack of foils is used the thickness of 
each foil should be small compared to the range of the recoil 
product to get good differential range curves. The concentration 
of the recoils in each layer or foil is determined suitably and, 
from a plot of concentration per layer versus the thickness 
penetrated, the range of the recoils can be obtained. 

The possible sources of error lie in the finite thickness 
of the target, the finite size of the collimator holes and the non- 
uniformity in the thickness of the foils used or in the thickness 

t 

of layers removed from the thick catcher. 


(ii) Thin target-thin catcher integral method - A thin target 
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(source) which is emitting recoils isotropically (like a fissile 

(source) is kept close to a stack of catcher foils. The activity 

At present beyond a thickness t, from the catcher surface facing 

the source is related to the range R, is given by Segre and 
27 

Wiegand/ by the equation 

= k (R-t) , 

where k is a constant. A plot of A^ versus t gives a straight 
line and the intercept on the t-axis gives a range R. Finkle 
lised this method for determining the ranges of 

235 

from U fission in alximiniiom catcher. 

The sources of error in this method are the finite thick- 
ness of the target (source) and the non-uni foimity in the thick- 
ness of the catcher foils. 

4.2.2 Thick Target Methods 

One of the main disadvantages of the thin target method is 

that they are useful only if the target material has a large 

nuclear reaction cross-section. If the target material has a 

low reaction-cross section, long periods of irradiation will be 

required. This difficulty can be overcome by using the thick- 

29 

target method developed by Sugarman et al . 

(i) Thick target-thick catcher method- A tajrget (source) , whose 
thickness is much larger than the range of the recoil product in 
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it# is sandwiched between two catcher foils. The catcher foils 
should be thick enough to stop all the recoils of interest. If 
the recoils are assumed to be emitted isotropically# a reasonable 
assumption if the target is amorphous and the projectiles are 
incident on it isotropically, the range R, of the recoil product 
in the target material is given by 



where is the activity of the recoil in one of the catchers# 

is the total activity of the recoil product both in the targe 

and the catcher foils# and t is the thickness of the tairget. 
on 

Niday has used this method for measuring the ranges of fission 

235 

products from the thermal neutron induced fission of U. 

31 

(ii) Thick target-thin catcher method - Alexander and Gazdik 

used targets that were not of negligible thickness but were able 

to obtain ranges of fission products in the catcher material. 

The target is sandwiched between stacks of catcher foils of 

different thickness. The fraction of the total yield of a 

particular kind of recoil product penetrating beyond a given 

^ 31 

thickness t of the external catcher is given by 



the target thickness and r is the average reciprocal 

R 


where W is 
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range in "the catcher material, c is a constant which relates 
the rate of velocity loss in the target material to the rate of 
velocity loss in the catcher material, and is given by the equa- 
tion; 


c' F 


t 


d W 



This method assumes that the recoil products are emitted isotro- 
pically. 


4.3 NUCLEAR EMULSIONS 

32 

Nuclear emulsions contain finely divided grains of silver 
bromide dispersed uniformly in gelatin. They differ from ordi- 
nary photographic emulsions in the following important aspects: 

i) Nuclear emulsions contain a high concentration of silver 
bromide, as much as four times that in ordinary photographic ernul- 
sions. 

ii) In nuclear emulsions the silver bromide grains are 
smaller in size and well separated compared to those in photo- 
graphic emulsion. Hence, unlike photogr aphic emulsions, nuclear 
emulsions can be used for registering individual events. 

iii) Nuclear emulsions are quite insensitive to light. 

The composition of the elements silver, brcxnine, iodine, 
hydrogen, carbon, nitrogen, oxygen and sulfur present in nuclear 
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emulsions can vary depending upon the make and the specific use. 

In addition to these, it is possible to load nuclear emulsions 
with other elements like lithium, boron, bismuth and uranitam for 
specific studies. 

An ionising particle passing through nuclear emulsion 
leaves a track containing a number of developable silver bromide 
grains. Emulsions with the smallest grain size are sensitive only 
to densely ionising particles like fission fragments while those 
with coarse grains can register tracks of <=( -particles, protons, 
mesons and electrons also. Thus, a proper choice of grain size 
allows a discrimination between various particles. 

The ranges or path lengths of ionising particles in nuclear 
emulsion, after development, can be measured using an optical 
microscope. The development of nuclear emulsions differ consider- 
ably, in many respects, from the techniques used for developing 
photographic emulsions mainly because nuclear emulsions plates 
are very thick and contain large quantities of silver. For thick 
emulsion plates, which are used particularly in high energy 

33 

studies, the "temperature method" of Dilworth ^ is 

particularly useful for obtaining uniform development of the 
whole emulsion without overdeveloping the outer layer. In this 
method, nuclear emulsion plates are first soaked in the developer 
solution for a long time at a terrperature low enough to prevent 
the development process. Subsequently, the temperature is raised , 
whereupon development proceeds uniformly throughout the emulsion. | 
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In using nuclear emulsions for range measurement, special 
precautions have to be taken to see that there is not even the 
slightest distortion of gelatin due to the development process. 
Also, under extreme experimental conditions, the composition of 
the emulsions may change and proper corrections have to be 
applied to the measured ranges. 

The energies of the ionising particles can be obtained 
from their ranges in nuclear emulsions using suitable range- 
energy relations. Though a number of range-energy relations 
-have been developed, they lack general applicability, and in many 
cases require arbitrary corrections, A part of the present work 
is concerned with the development of general range— energy rela- 
tions for heavy ions in complex media, including nuclear emul- 
sions and is discussed in detail in Chapter 3. 


4.4 SECTIONING TECHNIQUES 

The measurement of ranges of heavy ions like fission 
products, which are of the order of a fev/ microns in a given 
catcher, is not an easy task. To obtain proper range distri- 
butions sufficient data must be available and it is necessary to 
know concentrations per layer in depth increments of fractions 
of a micron. Though a large number of techniques have been deve- 
loped for making such measurements the ones involving the removal 
of thin layers from the surface of the catcher coupled with the 
measurement of implanted at<xns removed or remaining have proved 
the most successful. 
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4.4.1 Chemical Dissolution 

In this method suitable chemical etchants are used for 
removing thin layers from the surface of the catcher by chemical 
dissolution. The method has been used in the study of the range 
profiles in sodi\am chloride and potassium bromide^"^ and in GaAs.^^ 
The preferential etching rate on different crystallographic faces 
of the single crystal, at grain boundaries, dislocations, and, 
infact, at any point where local structure differs from the 
matrix, as for example, the implanted site, is a major disadvan- 
tage of this method. These difficulties, however, may not arise 
if amorphous catchers are used, 

4.4.2 Anodic Oxidation / Chemical Dissolution 

36 

In 1960 Davies ^ ad. developed a novel peeling tech- 
nique which has been responsible for much of the progress that 
has been achieved during the succeeding years in the study of 
range profiles, particularly of low energy ions. In this n^thod 
the speciman is made the anode of an electrolytic cell and anodiz- 
ed to a known thickness of oxide layer. The thickness of the 
oxide layer is controlled by the voltage applied. The oxide 
layer is subsequently selectively dissolved using suitable sol- 
vents, Though undoubtedly the most sensitive and reproducible, 

37—39 36 40 

it is at present limited to silicon, aliuninium, tungsten, 

gQTd, molybdenum^^ and indium antimonide. The limitation is 
not so much due to the difficulty in anodizing but in finding 
proper chemicals which selectively dissolve the oxide layer 
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without attacking the metal. 

^ *^o^^o sive Film / Chemical Dissolution. 

This method is based on the observation that when soft 
metals are immersed in an organic solvent containing halogen, a 
corrosive film is formed on the surface of the metal. This film 
can be selectively dissolved by chemicals which form complexes 
with metal halides. Andersen and SjzJrensen"^^ have made range 
studies in copper, silver and gold using this method. 

4.4.4 Mechanical Polishing 

This technique has been used to measure range profiles in 

45 

a wide variety of materials. The main disadvantages of this 
method is that the deformations introduced in the specimen due 
to mechanical abrasion can result in anomalously shortened ranges. 
Nonetheless, the universal applicability compensates for this 
drawback. 

4.4.5 Low Energy Ion Sputtering 

The ejection of atoms from the surface of the material 
bombarded with energetic ions is known as sputtering. In princi- 
ple, this method is applicable to all materials, but in practice, 
is very sensitive to changes in crystal stiructure. Since the 
sputtering rates are different for different substances there is 
a possibility of the implanted atoms being sputtered out at a 
faster or slower rate than the matrix. Of more concern in range 
profile studies is that in materials damaged by ion in^lantation. 
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where the damage profile changes along the path, a continuous, 
hut nonuniform change in sputtering rates through the range 
profile is likely. Not withstanding these drawbacks, the tech- 

. 47 

nique has been used by Lutz and Sizmann, and by P&hlaue et al . 
for the measurement of range profiles in copper and gallium 
arsenide, respectively, 

4.4.6 Electrolytic Dissolution 

It has been well known that metals can be quantitatively 
dissolved as the anode of an electrolytic cell. However, there 
have been no attempts to use this technique for range measurement. 
In the present work, differential ranges of the fission products 
^^Zr and ®^Sr (from fission) in copper have been 

measured using this technique. The experimental procedure is 
discussed in detail in Chapter 5. 



CHAPTER 5 


COULOMETRIC METHOD OF 
RANGE MEASUREMENT 


5 .1 INTRODUCTION 

Fission products from the thermal neutron induced fission 

235 

of U/ after collimation, were allowed to penetrate a copper 
foil at ri^t angles to its surface. Thin layers of accurately 
known thicknesses were then electrolytically dissolved out from 
the catcher foil. The thickness of the layers could be controll 
ed by varying the duration of electrolysis and the current 
strength. The fission products ^ose range was to be measured 
was separated by radiochemical methods frcxn the electrolyte 
after the dissolution of each layer. The range was obtained 
from a plot cf the activity of the radioactive fission product 
in each layer acainst the depth at which the layer is situated • 
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5 . 2 TARGET PREPARATION 

For accurate differential range itieasurements ’thin' targets 
required so that there is negligible energy— loss of the 
fission product in the target material itself. 

Thin targets of uranium can he obtained by its electrodepo" 
sition as hydrous oxide at the cathode. The electrolytic cell 
is shown in Fig. 5,1. it consists of a glass cylinder G,'^3 cm 
in diameter andf--^6 cm in height, with a flange at the lower end. 

The copper foil C,~(2 cm x 2 cm) in area, on ^ich uranium hydro- 
xide is to be deposited is placed on a copper disc D. A copper 
wire is soldered to D for making electrical contact. The area 
of the deposit on the copper foil C v^ich acts as the cathode 
was controlled by meanSof a teflon disc T, having a punched out 
hole of 1.5 cm diameter. The flange of G rests on T and is held 
in position tightly with the help of metal rings F, F' and screws 
v/hich rnake the lovrer part leak-proof. The copper disc is insitlated 
from the lower metal ring by mica. The platiniim disc P, <-'2 cm 
in diameter, acts as the anode. 

The electrolyte consisted of equal volimes of 0.2 M perchlo- 
ric acid and 0,15 M ammonium formate. To about 20 ml of the 
electrolyte in the electrolytic cell 0,5 ml of a dilute nitric 
acid solution of uranyl nitrate (containing — 200 Pg of uranium 
enriched in (^^60%)) was added. The anode and the cathode 

were about 2 cm apart. Electrolysis was carried out at room 
temperature (<^25^0) . A current density of "—SO mil li amperes 
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per square centimeter of the cathode vras maintained during 
electro lys is . 

Under these experimental conditions anodic oxidation of 

4R 

the formate ion proceeds as follows: 

HCOO" ^ CO 2 + + 2 e" . 

at E° = 0.196 volt. 

This reaction produces half a hydrogen ion equivalent at the 
anode while the cathode reaction, involving hydrogen ion reduc- 
tion, consxmes one hydrogen ion equivalent per faraday of elec- 
tricity. The net cell reaction causes the loss of half a hydrogen 
ion equivalent per faraday of electricity. The acidity of the 
Electrolyte is slowly reduced and tlnis causes t±ie quantitative 
deposition of uranium hydroxide at the cathode because of the 
change in the pH. Electrolysis was carried out till approxima- 

23 ^ 

tely 75 fj>g of ''U v/as deposited. After this the cathode with 
the uranium hydroxide deposit vras removed, washed with distilled 
water and dried at 110*^0. 

5.3 TARGET ARRANGEMENT AND IRRADIATION 

The deposited uranium on copper backing was then incorpo- 
rated into the target assembly for the purpose of thermal -neutron 
irradiation. Figure 5.2 shows the target assenbly. The uranium 
target on copper backing T, and the copper catcher foil C were 
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separated by the collimator M, The colliinator was a 5 mm-thick 

aluminium block with about 100 parallel holes of 0,8 mm diameter 

in an area of '^2.5 square centimeter. This assembly was then 

wrapped in a pure aluminium foil, sealed in a glass container 

evacuated to a pressure of'-'5 x 10 meters of mercury'', enclosed 

in water tight alxjminium can and irradiated in the Apsara reactor, 

Bhabha Atomic Research Centre, Bombay, with a thermal neutron 

13 2 

flxix of the order of 10 neutrons/cm sec. To get sufficient 
activity for good counting statistics irradiation was done for 
seven days after vAiich the sanples were allowed to "cool" for a 
few days so that the short-lived activation products decayed 
away and the samples could be handled v/ith safety. The assembly 
was then opened and the catcher foils were subjected to the treat- 
ment described in the succeeding sections. 

5 .4 REMOVAL OF THIN LAYERS FROM THE CATCHER FOIL 

Since copper has a low oxidation potential, thin layers 
could be removed from the copper catcher in an electrolytic cell. 
From the total quantity of electricity passed during electrolysis 
and the area of the anode in contact with the electsrolyte the 
thickness of the layer dissolved can be calculated. If there is any 
purely chemical dissolution of the anode then corrections would 
have to be applied to the thickness calculated on the basis of 
purely electrochemical dissolution. 
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5.4.1 Electrolytic Cell 

The electrolytic cell is shown in Fig, 5,3, It consists of 
a rectangular perspex chamber of dimensions 2 x3 >:5 cm. The copper 
catcher foil from V7hich thin layers are to be dissolved acts as 
the anode . The cathode consists of two identical copper 
strips having the same dimensions as the anode and mounted syirme- 
trically on either side of the anode at a distance of 1,5 cm 
from each other. The smaller the distance between the anode and 
the cathode compared to their own dimensions, the more uniform 
would be the electrical field between them resulting in more ixni- 
form dissolution of the anode surface. The anode and the cathode 
foils are fixed rigidly on a perspex frame so that during 
the removal of successive layers the geometry of the electrodes 
could be kept unaltered. 

A Sargent coulometric current source, model-IV (supplied 
by E.H. Sargent and Co., U.S.A.) was used to obtain a constant 
current of precisely known strength for electrolysis. Six pre- 
determined constant currents in the range of 4.825 millianpers to 

193.0 milliampers could be selected. The timer v^^ich measures 
1 til 

time in -j-— of a second is automatically switched on when the 
cell current is applied. The cell resistance of— 250 ohms was in 
the specified range for which current regulation was +0.1 per cent 

5.4.2 Choice of Electrolyte 

A 10 % solution of sodiiam nitrate in 0,5 N nitric acid was 
used as the electrolyte. The choice of this electrolyte was 
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Fig. 5.3 Electrolytic Cell 
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based on the following considerations: 

A. Current Efficiency ; The current efficiency of the electrolyte 
is one of the most important aspects that has to be taken into 
consideration in any constant current coulometric experiment. The 
current efficiency of the above electrolyte was studied before- 
hand as described below; 

In the electrolytic cell a platinum wire was used as the 
cathode, the other electrode being a thin copper strip. Electro- 
lysis was carried out for a short duration ('^3 seconds) using a 
constant current of 4.825 milliampers. From the quantity of the 
electricity passed the quantity of copper expected to be deposi- 
ted on the cathode was calculated on the basis of Faraday* s law 
of electrolysis. The actual amount of copper deposited was deter- 

49 

mined spectrophotometrically by the dithizone method as follows: 

Calibration curve - A 0,01?^ copper solution was prepared 
by dissolvinc 0,1969 g of ijneff loresced crystals of copper sul- 
fate (Cu 50^.5H20) in 500 ml of 0,1 N sulfuric acid. Frcxn this 
solution weaker standards could be prepared by dilution with 0.1 N 
sulfuric acid. 10 ml of standard solution was taken in a small 
(50 ml) separatory funnel. 0.001% dithizone solution in carbon 
tetrachloride was delivered in small quantities from a micro- 
burette, shaking well the contents of the separatory funnel after 
each addition. The quantity of dithizone used was sufficient 
(i.e., slightly in excess) to extract the copper completely as 
copper dithizonate into carbon tetrachloride. This was indicated 
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by the change in the color of the final solution from that of 
pure copper dithizonate (red-violet). The exact voliime of 
dithizone solution used was noted. 

A few drops of carbon tetrachloride was allowed to flow 

out of the fmnel to displace any aqueous solution in the bore 

of the stop-cock and the stem was dried with filter paper rolled 

around a glass rod. The copper dithizonate solution was then 

delivered into the quartz cell and the per cent transmittancy, T, 

of the solution was determined at 510 nm using spectronic-20 

colorimeter (Bauseh & Lorrib Optical Co,, U. The procedure 

was repeated with different standard solutions of copper. The 

calibration curve is shown in Pia. 5.4 in which loa (1/T) is 

'10 

plotted against the concentration of copper ( /ug/l ) . 

f 

Determination of the Weight of Copper Deposited on the Cathode -■ 
The copper deposited on the platintim cathode was dissolved in 
dilute nitric acid (2N HNO^) . The solution was evaporated to 
dryness and heated strongly to remove the oxides of nit3X>gen 
completely. The residue was dissolved in 5 mi of 0.1 F sulfuric 
acid. The 'vdaole of this solution was then shaken with 0.001% 
dithizone solution in a 50 ml separatory funnel as described 
before to extract copper completely into carbon tetrachloride as 
copper dithizonate. The exact volume of dithizone solution used 
was noted. The copper dithizonate solution was suitably diluted 
and the transmittancy was determined at 510 nm. The concentra- 
tion of copper solution was obtained from the calibration curve 
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and from this the wei^t of copper deposit was calculated. From 
the values given in Table 5.1 it can be seen that the current 
efficiency of the electrolyte can be taken to be 100% without 
introducing significant errors in range measurement. 

B. Chemical Dissolution of Copper : The coulometric current 
source gives the total quantity of electricity passed during 
electrolysis and hence the amount of the anode material that has 
been dissolved electrolytically. if there is any purely chemical 
dissolution it is essential to account for it in calculating the 
thickness of the layers dissolved out. The following experiment 
was performed for this. 

The copper anode was weighed accurately and the electro- 
lysis was carried out for a known duration of time, approxima- 
tely one hour. The loss in weight of the anode after electro- 
Ivsis cave the weight of copper that had gone into solution both 
due to electrolysis and any possible chemical etching. As the 
current efficiency of the elect3x>lyte had been fotand to be lOO 
per cent^ the v\reight of copper dissolved electrolytically could 
be calculated from the quantity of electricity passed. Any 
difference between these two weights should give the weight of 
copper that had gone into solution due to purely chemical disso- 
lution^ and knowing the duration of electrolysis the rate of 
chemical attack could be calculated. Table 5.2 gives the experi- 
mental results. Since the difference in weights is within the 
limit of weighing accuracy (all weighings were done using a 



letween the weights of copper deposit determined by the dithizone 
and that calculated using Faradays law of electrolysis (M^) 
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T^le 5.2 


Comparison between the wei^ts of copper dissolved from the 
anode as determined by weighing and that calculated using 

Faraday's law of electrolysis 


Total quan- 

Wt. of the 

Wt. of the 

Wt, of the 

Loss in wt. 

tity of ele- 

anode mate- 

anode before 

anode after 

of the anode 

ctricity 

rial expec- 

electrolysis 

elect roly- 

material 

passed (MP ) 

ted to dis- 
solve (mg) 

( g) 

sis { g ) 

(mg) 

231.65 

7.36 

2.91238 

2.90501 

7.37 

313.10 

9.95 

2.90501 

2,89491 

10.10 


balance which can weigh down to 0,01 mg) no corrections for chemi- 
cal attack was deemed necessary in calculating the thickness of 
the layers dissolved. 


C. Ease of Chemical Separation of Fission Products ; since the 
electrolyte contains only sodium nitrate and nitric acid separa- 
tion of fission products from it after electrolysis could be done 

easily using well tested methods available in literature with 

50, 51, 52 

only minor modifications if, at all, required. 


5,4.3 Anodic Removal of Thin Lavers from the Copper Catcher 

The electrolytic cell was set up as described in Sec. 5. 4,1 

with copper catcher foil (5x1.5 cm) as tlie anode. 15 ml of the 
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electrolyte was added to the cell covering the anode to a height 
of approximately 2,5 cm. Electrolysis was carried out for an 
accurately known duration of the order of lOO seconds. The dura- 
tion of the electrolysis was chosen by taking into consideration 
the theoretical range of the fission product so that practically 
all the fission product activities could be dissolved out in 
about 35 layers of catcher in order to obtain a good differential 
range curve. The current was maintained constant at 48,25 milli- 

amperes corresponding to a current density of approximately 10 
2 

mini amperes/cm of the anode area. After each electrolysis the 
electrodes were carefully washed with demineralized water vhichi 
was collected inside the cell. The solution in the electrolytic 
cell containing the fission products from the dissolved layer was 
transferred to a labelled 50 ml centrifuge tvibe. The electrolytic 
cell was washed with demineralized water and the washings were 
added to the same centrifuge tube. Electrolysis was repeated 
35-40 times as described above and the solutions containing the 
fission product from each layer were stored separately in labelled 
50 ml centrifuge tubes. The fission products vhose ranges are 
to be determined were then separated from the above solutions 
radiochemiically as described in the next section. 

5.5 RADIOCHEMICAL SEPARATION OF FISSION PRODUCTS MJD COTOTTIHG 

Radiochemical separations were done by "carrier” procedures 
in \diich milligram amounts of the elements of interest are added 
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at the beginning of the procedure. Since the amount of inactive 
carrier added is known, a determination of the amount present at 
the end of separation gives the chemical yield of the seoaration 
process and an appropriate correction is applied to the observed 
radioactivity to obtain the actual activity. To achieve aood 
radiochemical purity one has to resort to repeated purification 
involving several precipitation and scavenging operations. Also,, 
the procedures are often quite time-consimiing. Talcing these 
factors into consideration, it was advantageous to choose fission 
products with reasonably long half life and hi^ fission yield. 
The fission products chosen for study in this vork are given in 
Table 5.3 along with their half-life and cumulative yield which 
represents the probability of formation in fission of a nuclide 
of given mass nunber after prompt and delayed neutron enission 
and the decay of its short-lived precursor nuclides. 


Table 5.3 

Half -life and the cxmulative yields for fission 
products chosen for range measurenoats 


Fission 

product 

Half-life 
( days) 

Cumulative yield 
{%) 

140„ 

12.8 



6.25 

40 

65.0 

6.53 

®"sr 

38 

51.0 

4.73 
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5,5,1 Bari\iin-140 

Baritim can be separated frcm the other fission products 

by the specific precipitation at low temperatiire as Baci 2 . 2 H 20 

by means of concentrated hydrochloric acid- ethyl ether mixture. 

The procedure outlined below, which is similar to that developed 
50 

by Glendenin, consists of three precipitations as chloride, 
the first and second being followed by ferric hydroxide scaven- 
ging purification steps. The final precipitation of Barium as 
Baritun chromate (BaCrO^) was done in the presence of inactive 
strontixam as hold-back carrier. The different steps involved in 
the separation procedvire are described below: 

i) 2 ml of barium chloride solution containing 20 mg of 

2 + 

Ba was added to the whole of the electrolyte solution contain- 
ing the fission products from the layer dissolved from the 
catcher. 

ii) Barixom was precipitated as barixam carbonate by adding 
saturated sodixm carbonate solution in slight excess and centri- 
fuged. The supemate was discarded. 

iii) The barixam carbonate precipitate was dissolved in 30 mJ 
of ice-cold etheir-hydrochloric acid mixture containing 1 volxame 
diethyl ether and 5 volxames of concentrated hydrochloric acid, 
and cooled in ice with occasional stirring. 'White precipitate of 
barixam chloride separates out. The precipitate was centrifuged 
and washed thrice with 5 ml poartions of ice-cold ether-hydro- 
chloric acid mixture. 


4 
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iv) The baritun chloride precipitate was dissolved in 5 ml 

3 + 

of deionised water, 5 ml of Fe scavenger, prepared by dissol- 
ving ferric c±iloride (FeCl 2 . 8 H 20 ) in very dilute hydrochloric 
acid was added. Ferric hydroxide v;as precipitated completely 
by adding concentrated ammonia dropwise with constant stirring 
and centrifuged. 

v) Steps (ii), (iii) and (iv) were repeated with the super- 
nate from step (iv) . 

vi) From the supemate in step (v) barium chloride vras 
precipitated via steps (ii) and (iii) . 

vii) The barium chloride precipitate was dissolved in 5 ml 
of water. 20 mg of Sr^"^ carrier, prepared by dissolving strontium 
chloride (SrCl 2 - 2 H 20 ) in water, was added to act as hold bach 
carrier for strontium isotopes. 

viii) The solution was neutralized with anmonia till just 
alhaline. 6 n hydrochloric acid was added dropwise to make the 
solution just acidic to methyl orange. To this solution 1 ml of 
6 M acetic acid and 2 ml of 6 M ammonium acetate were added- The 
contents were heated to nearly boiling and 1 ml of 1.5 M lootassium 
chromate was added dropwise with stirring to precipitate barium 
chromate. Heating v;as continued for two minutes. 

ix) A f/Jhatman No. 42 filter paper disc, '^2.5 cm in diameter, 
was washed with distilled water and ethyl alcohol, dried at llo C 
for 15 minutes, cooled to room teirperature and weighed. 
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x) The barivim chromate (BaCro_^) precipitate was filtered 
on-to above filter paper disc in a Hirsch chinney. The precipi- 
tate was washed thrice with 5 ml portions of distilled water/ 
thrice with 5 ml portions of ethyl alcohol, dried at llo'^C for 

15 minutes, cooled to room temperature and weighed. The chemical 
yield of barium chromate was '~70%. 

xi) The barium chromate precipitate was mounted on an 
aluminixim holder and the radioactivity was measured using a 
sharp low background beta proportional counter. 

5,5.3 Zirconium-95 

The separation of zirconium was done b^'^ using the procedure 

51 

developed by Hahn and Skonieczny. In this method, after the 
initial decontamination by scavenging with lanthanum fluoride 
in the presence of inactive zirconium (Zr ) carrier, zirconium 
is separated as insoluble barium fluozirc.nate. The purification 
of zirconium by two more precipitations as barium -fltio zircon ate 
in acid mediiim is followed by its conversion to zirconium tetra— 
mandelate in vAiich form it is weighed and counted. The separa- 
tion procedure is given below: 

i) To the electrolyte containing the fission products 2 ml 

4 + 

of zirconium nitrate solution containing 5 mg of Zr was added. 
The contents were transferred to a lusteroid tube. To this 0.5 ml 
27 N hydrofluoric acid was added and mixed thoroughly. 

ii) To the above solution 0.25 ml of ' lanthanum nitrate 
solution containing 10 mg of La /ml was added, stirred well. 
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and. * cent ri fug ed briefly. A second 0.25 ml of lanthanum nitrate 
solution was added and the precipitate v;ras centrifuged down on 
top of the first precipitate. The supemate was decanted into 
another lusteroid tiibe. 

iii) To the above supemate 1 ml of Ba ( 1302)2 solution contain- 
. 2 + 

ing 50 mg of Ba per ml was added, stirred well, and allowed to 
stand for a minute. The barium fluozirconate precipitate was centri- 
fuged and the supemate discarded. 

iv) To the lusteroid tube containing barium fluozirconate 
precipitate 2 ml of distilled water and 1 ml of saturated boric 
acid solution were added and stirred well. The precipitate was 
dissolved by adding 0,25 ml of concentrated nitric acid and heat- 
ing on a waterbath. To the clear solution 0,5 ml of barixam 
nitrate solution was added and mixed vjell. 0.25 ml of concentra- 
ted hydrofluoric acid was added to precipitate zirconium as barium 
fluozirconate. The precipitate was centrifuged and the supemate 
discarded, 

v) Step ( iv) was repeated with the precipitate. 

vi) The precipitate from step (v) was suspended in 2 ml of 
water. 0,5 ml of saturated boric acid solution and 2 ml of 6w 
hydrochloric acid was added and stirred well. The solution was 
made distinctly alkaline by adding 6N sodixom hydroxide x-xith 
constant stirring. The zlrconixxm hydroxide precipitate v/as centri- 
fuged and the supemate was discarded. 

vii) The zirconium hydroxide precipitate was dissolved in 
3 ml of concentrated nitric acid and diluted to 6 ml. The 
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solution was transferred to a glass centrifuge tube. The luste— 
roid tube was rinsed twice with 1 ml portions of water and the 
washings were combined with the solution. 10 ml of saturated 
mandelic acid solution in water was added, stirred well and heated 
on a waterbath at 80°-90°C for about 20 minutes. Zirconixim tetra- 
mandelate Zr( CgH^CHOHCOO) ^ precipitated out. 

viii) The contents were filtered hot on-to a ^#iatrnan wq. 42 
filter paper disc, '^ 2,5 cm in diameter, -vihich had been previously 
washed with water, ethyl alcohol and diethyl ether, dried at 50®C 
for 15 minutes cooled to room tenperature and weighed. The pre- 
cipitate was washed with 10 ml of 5% mandelic acid in 2 % hydro- 
chloric acid followed by three 5 ml portions of diethyl ether. 

The precipitate was dried for 15 minutes at 50*^C, cooled to room 
temperature 3nd weighed. The chemical yield of zirconium tetra— 
mandelate was^40%. 

ix) The zirconiiom tetramandelate precipitate was then 

counted using a "sharp' low background beta proportional counter. 

Counting was done immediately after separation and no corrections 

95 

have been applied to the contributions frcm the daughter SB 
which has a half-life of 35 days, 

5,5.3 Strontium-89 

The procedure for the separation of strontium from fission 

5 2 . • 

products was developed by Beaufait and Lukens. Strontium is 
first precipitated as nitrate at low temperature by means of 
yellow fuming nitric acid. This is an excellent decontamination 
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step, the major impurity being barium. The strontium nitrate 
precipitate is purified by ferric hydroxide scavenging followed 
by quanitative precipitation of baritun as chromate in the pre- 
sence of inactive baritm carrier, .^-t this stage an oxalate 
precipitation followed by a nitrate precipitation will give the 
necessary radiochemical purity. The details of the method are 
given below: 

i) To the sample solution containing the fission products 

2 + 

2 ml of strontium chloride solution containing 20 mg of Sr was 
added. The solution was concoitrated to a volume of about 2 ml. 

ii) 20 ml of fuming nitric acid was added to the above 
solution and cooled in an ice-bath for 20 minutes with occassional 
stirring. The strontium nitrate precipitate was settled by centri- 
fuging at low speed for about 10 minutes and the supemate was 
decanted off. 

iii) The stronti\am nitrate precipitation was dissolved in 

3 + 

7 ml of distilled water 5 mg of Fe solution was added to the 
above solution. Ferric hydroxide was precipitated completely by 
adding 6N ammoniim hydroxide drcpwise with stirring. The ferric 
hydroxide scavenge was centrifuged and the char solution was 
transferred to another clean centrifuge tube. The ferric hydrox- 
ide scanvenge was washed by slurrying with 7 ml of 5°^ ammonium 
nitrate solution containing 2—3 drops of 6 N ammonium hydroxide 
and again centrifuged. The clear washings were combined with the 
above solution. 
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2*f 

iv) 20 mg of Ba carrier, prepared by dissolving barium 
chloride (Baci 2 . 2 H 20 ) in water, was added to the above solution. 

The solution was neutralized with 6N nitric acid, checking the 
pH with a pH papen 1 ml of 6 M acetic acid and 2 ml of 6 M ansno- 
nium acetate were then added to the solution and it was heated 
nearly to boiling. Barium chromate was completely precipitated 
from this solution by adding 1 ml of 1.5 M potass ixim chromate 
dropwise with constant stirring. After centrifuging the clear 
supemate was transferred to another centrifuge tube. 

To the above solution 2 ml of 6N ammonium hydroxide was 
added and heated nearly to boiling. Strontiijm was precipitated 
as oxalate by adding 5 ml of saturated ammonium oxalate solution 
and centrifuged. 

vi) The strontium oxalate precipitate was dissolved in 4 ml 
of 6N nitric acid. 15 ml of yellow fuming nitric acid was added 
to it and cooled in ice for about 20 minutes with occasional stirr- 
ing, and then centrifuged. The supemate was discarded. 

vii) The strontium nitrate precipitate was dissolved in 
10 ml of distilled water, neutralized x^ith ammonium hydroxide 
and heated nearly to boiling. Strontium was precipitated as 
carbonate by adding 2 ml of 1 M sodium carbonate solution. Heating 
was continued for several minutes foliov/ed by cooling in ice for 
about 10 minutes. 

viii) A Wnatman No. 42 filter paper disc, <^2.5 cm diameter, 
was washed with distilled water and ethyl alcohol, dried at 110°C 
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for 15 minutes, cooled to room temperature and weighed. 

ix) The strontium carbonate precipitate was filtered on-to 
the above filter paper disc in a Hirsch chirmey. The precipi- 
tate was washed thrice with 5 ml portions of water and thrice 
with 5 ml portions of ethyl alcohol. The precipitate was dried 
at 110°C for 10 minutes, cooled to room tOTperature and weighed. 
The chemical yield of strontium carbonate was 60%. 

x) The strontium carbonate precipitate was mounted on a 
alxaminium holder and the radioactivity was measured in a sharp 
low background beta proportional counter. 

5.6 RESULTS 

After the counting the activity in the different layers 

the count rates were normalized to the same weirtt of the carrier 

and the same arbitrairy initial time. The thickness of the layers 

dissolved out from the copper catcher and the normalized activity 

in them are shown in Tables 5.4, 5.5 and 5.6. Figures 5.5, 5.6 

and 5.7 show the plot of the normalized count rate of different 

layers against the depth of the mid-point of the layer from the 

exposed surface. The well resolved, almost gaussian shaped, 

peaks allow the most probable ranges to be ascertained with good 

140 95 , 89 

accuracy. The most probable ranges of Ba, Zr and Sr are 
shown along with the full width at half maximum for the curves 


in Table 5.7. 
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Table 5«4 

Thickness of the layers dissolved out from the copper catcher 

and the relative radioactivity (count rate normalized to 100% 

chemical yield and the same arbitrary initial tin®) of 

in them, A constant current of 48,25 milliangjeres/seccMid was 

used for electrolysis. The total area (both faces) of the 

copper catcher (anode) in contact with the electrolyte was 
2 

9,96 cm , 


Layer No, 

Duration of 
electrolysis 
(seconds) 

Thickness of 
the layer 
dissolved 

( pg/cn?) 

Relative 

radioactivity 

(counts/ndnutes) 

1 

5 

3 

4 

1* 

100.5 

178.2 


2 

100.0 

177.3 

46.21 + 1.29 

3 

100.1 

177.5 

56.56 + 1.70 

4* 

100.0 

177.3 


5 

100.0 

177.3 

47.13 + 1.26 

6 

100.0 

177.3 

27.13 + 0.97 

7 

100.0 

177.3 

48.70 + 1.43 

8 

99,9 

177.1 


9 

100.0 

177.3 

45.20 + 1.38 

10 

100.0 

177.3 


11 

99.9 

177.1 

38.32 + 1.03 

12 

100.1 

177.5 

31.99 Hr 0.99 

13 

100.0 

177.3 

56.52 + 1.69 

14 

100.4 

178,1 

65.56 + 2.05 

15 

100.0 

177.3 

69.15 + 2.02 


contd 


Table 5.4 (contd.) 
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1 

2 

3 

4 

16 

100,1 

177.5 

95.39 + 2.99 

17. 

99.9 

177.1 

108.29 + 3.37 

18. 

100.0 

177.3 

128.38 + 3.81 

19 

99.9 

177.1 

140.71 + 4.34 

20 

100.0 

177.3 

158.22 + 4.91 

21. 

100,2 

177.6 

171.57 + 4.96 

22 

100.1 

177,5 

185.36 + 5.59 

23 

100.1 

177.5 

187.65 + 5.69 

24 

100.1 

177.5 

181.09 + 5.51 

25 

100.2 

177.6 

165.96 + 4.92 

26 

100.0 

177.3 

152.58 + 4.74 

27 

100.1 

177.5 

130.37 + 3.99 

28 

99.8 

176.9 

103.30 + 3.56 

29 

100.0 

177.3 

95.10 + 2.98 

30 

100.1 

177.5 

74.57 + 2.22 

31 

99.9 

177.1 

59.69 + 1.80 

32 

100,0 

177.3 

43.36 + 1.35 

33 

100.1 

177.5 

33.52 + 1.12 

34 

100.0 

177.3 

27,50 + 0.83 

35 

99.9 

177.1 

20.50 + 0,64 

36 

100.4 

178.1 

18.93 + 0.57 

37 

100.0 

177.3 

15.33 + 0.44 

38 

100.0 

177.3 

10.37 + 0.32 

39 

100.2 

177.6 

8.99 + 0.26 


*These samples were spoiled during chemical separation. 
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Table 5.5 


Thickness of the layers dissolved out from the copper catcher 

and the relative radioactivity (i.e., count rate normalized to 

100% chemical yield and the same arbitrary initial time) of 
95 

Zr in them, A constant current of 48.25 mamp/sec was used 
for electrolysis. The total area (both faces) of the copper 


catcher 

(anode) in contact 

2 

with the electrolyte was 6.08 cm , 

Layer No, 

Duration of 
electrolysi s 
(seconds) 

Thickness of 
the layer 
dissolved 

(fjLg/cm^) 

Relative 

radioactivity 

(counts/minutes) 

1 

2 

3 

4 

1 

150,0 

391.8 

39.73 + 1.21 

2 

150.2 

392.3 

22.02 + 0.66 

3 

130.1 

355.9 

20.32 + 0.62 

4 

119.7 

313.2 

26.10 + 0.78 

5* 

120.0 

313.7 


6 

120.0 

313.7 

15.66 + 0.47 

7* 

120.1 

313.7 


8 

120.3 

314.2 

16.33 + 0.48 

9 

119.7 

313.2 


10 

119.9 

313.7 

22.01 + 0^65 

11 

120.0 

313.7 

30.58 + 0.85 

12 

120.0 

313.7 

37.58 + 1.17 

13 

120.0 

313,7 

42.59 + 1.06 

14 

120, 3 

314.2 

44.77 + 1. 14 


contd, 
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Table 5.5 (contd.) 


1 

2 

3 

4 

15 

120.0 

313.7 

89.37 + 2.22 

16 

120.0 

313.7 

82.09 + 2.24 

17 

119.8 

312.9 

101.39 + 3.11 

18 

120.0 

313.7 

118.39 + 3.70 

19 

120.0 

313.7 

117.86 + 3.21 

20 

120.0 

313.7 

101.68 + 4.01 

21 

120.0 

313.7 

75.82 + 1.96 

22 

120.0 

313.7 

36.72 ± 0.84 

23 

120.0 

313.7 

32.57 + 0.91 

24 

120.0 

313.7 

14.96 + 0.46 

25 

119.7 

313.2 

8.97 + 0.27 

26 

120.0 

313.7 

12.14 + 0.38 

27* 

120.3 

314.20 


28 

120.1 

313.7 

2.16 + 0.22 

*These 

samples v/ere 

spoiled during chemical 

separation 
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Table 5«6 

Thickness of the layers dissolved out from the copper catcher 

and the relative radioactivity (count rate normalized to 100% 

89 

chemical yield and the same arbitrary initial time) of Sr 
in them. A constant current of 48,25 mamp/sec was used for 
electrolysis. The total area (both faces) of the copper 
catcher (anode) in contact with the electrolyte was 7.84 cm , 


Layer No, 

Duration of 
electrolysis 
(seconds) 

Thickness of 
the layer 
dissolved 

(/*g/cm^) 

Relative 
radioactivity 
( c ounts/minutes) 

1 

2 

3 

4 

1 

300.0 

607.9 

47.17 + 1.49 

2* 

300.0 

607.9 


3 

200.0 

405.2 

51.67 + 1.57 

4 

152.2 

308.4 

47.19 + 1.48 

5 

148.0 

299.5 

49.34 + 1.52 

6 

150.0 

303.9 

35.94 + 1.02 

7 

150.2 

304.3 

31.90 + 0.72 

8 

150.0 

303.9 

39.23 + 1.23 

9 

150.0 

303.9 

48.64 + 1.33 

10 

150.0 

303.9 

40.46 + 1.00 

11 

150.0 

303.9 

48.90 + 1.51 

12 

150.0 

303.9 

46.04 + 1.36 

13 

150.1 

304.1 

55.59 + 1.68 

14 

150.0 

303.9 

65.80 + 2.05 


contd. 
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Table 5.6 (contd. ) 


1 

2 

3 

4 

15 

150.0 

303.9 

96.40 + 2.97 

16 

150.0 

303.9 

130.70 + 4.01 

17 

150.0 

303.9 

166.24 + 4.95 

18 

150.0 

303.9 

183.35 + 5.38 

19 

150.1 

304.1 

214.80 + 6.45 

20 

150.0 

303.9 

190.41 + 5.94 

21 

150.0 

303.9 


22 

150.2 

304.3 

113.81 + 46.2 

23 

150.1 

304.1 

63.38 + 1.08 

24 

150.2 

304.3 

16.16 + 0.50 

25* 

149.8 

299.9 


26 

150.0 

303.9 

13.26 + 0.39 

27 

150.0 

303.9 

11.18 + 0.34 

28* 

150.0 

303.9 


29 

150.1 

304.1 

16.63 + 0.49 

30 

150.8 

305.6 

11.63 + 0.31 

31 

150.0 

303.9 

6.61 + 0.21 

■* These 

samples were 

spoiled during chemical 

separation. 
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20 
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Penetration depth (mg/cm^) 


Fig. 5.6 Differential range distribution of the fission product 
in Copper. 





The most probable ranges R (mg/cm ) and the full width at Ralf- 

P 2 

maxima (mg/ cm ) for the peaks 


I 
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5 .7 SOURCES OF ERROR IN THE MEASUREMENT 

5.7*1 Determinate Errors in the Radioch^lcal Procedures 

Some problems usually encotintered in radiochemical work 
like determination of counting efficiency for the radiation 
and the standardization of the carrier solution did not appear 
in this study because only relative activities were needed. 

The primary sources of error to be considered lay in weighing, 
counting, radioactive impurity and scattering of beta radiation 
in the samples. 

Ii 

The samples were weighed using a balance which can read 
upto O.Ol mg. Since the samples weighed 60-80 mg the accuracy 
in weighing was of the order of 0 , 15 %. The chemical procedures 
adopted have been shown to give adequate radioactive purity 
which was confirmed by careful monitoring of the decay of some 
of the samples. Since for the same fission product the chonoi- 
cal yield for the different samples did not differ appreciably 
no corrections were rrquired for self absorption and back- 
scattering of the beta radiation. 

5»7*2 Non-Uniformitv in the Thickness of the Layers Dissolved 
from the Catcher Foil 

The electric field between the parallel electrodes is not 
uniform, the lines of force being more crowded near the edges. 
Hence the rate of dissolution of the anode near the edges is 
higher than at the centre of the anode. The error in range 
measurement due to this effect was minimized by using electrodes 
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with large area and keeping the distance between them small to 
ensure a more xmiform electric field between than. 

5.7.3 Activation Impurities 

140 95 89 

Production of Ba^. Zr and Sr by neut 2 X)n activation 

of impurities present in the catcher would cause considerable 

error in the range measurement of these fission products. To 

ensure that the nuclides whose ranges were measured were produc- 

235 

ed only by the fission of U a 'blank run' was performed- 

Irradiations were carried out as described earlier but without 
235 

the U target. After irradiation the catcher foil was analysed 

for the activities present. The gamma ray spectra did not show 

140 95 89 

any photopeaks for Ba, Zr and Sr. 

5.7.4 Collimation Angle 

Because of the finite diir.ensions of the collimator holes, 

the fission products can penetrate the catcher surface at angles 

other than 90*^. In the present work a 5 mm thick collimator 

having 0.8 mm diameter holes was used. The maximum deviation 

0 from 90° at which the fission product could penetrate the 

catcher surface, calculated using simple geometric considerations, 

was'^9°. If all the fission products were to penetrate the 

catcher surface at (90° +9 ) the measured ranaes will be 

— max 

shorter by ^ 1 . 2 % from the true ranges. However, in actual 
measurements, the error would be smaller than this since all 
penetration angles between 90° and (90° + equally 
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pjTobsble. The rsnaes oiven in Table 5.7 have not been coirirecte<3 
for this effect, 

5.8 DISCUSSION 

5.8-.1 Comparison with Other Methods 

No data on ranges of the fission products ^^Zr and 

89 

Sr in popper are available in literature. Hence, a comparison 
of the presQit measurement with other methods used for range 
measurements has not been possible. Ws could only cxjmpare our 
experimental values with the corresponding theoretical values 
calculated by means of a reliable range-energy equation v^ich 
has beei tested in the case of other media. This is given in 
the next section. 


5.8.2 Comparison v;ith Calculated Ranges 

The ranges of the fission products can be calculated using 
4 

the equation: 


R = 




127.3 f(Z2) t4.7622 + f{Z^)] 


(5.1) 


where R is the range in mg/cm^ of the fission product of mass 
and atomic nuinbers and respectively; A 2 Z 2 are the 
mass and atomic numbers of the catcher material (stopping 

medium) ; V. is the initial velocity of the fission product and 

9 8 

V = e'^/fi ; both and are repressed in \inits of 10 cm/sec. 

4 

The f(Z) values are given by the Equation: 
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f(Z) = 0.28 for Z 4 : 45.5 , (5.2) 

f(Z) = Z^/^ for Z > 45.5 . (5.3) 

These f(z) values were derived using the experimental ranges of 
fission products. 

(a) Fission product velocities - For calculating the ranges of 

the fission products their initial velocities are required- 

From the experimentally measured kinetic energi'=>s of the primary 

fission fragments the velocities of the fission products were 

calculated as follows. For the given fission- product of mass 

number , the precursor fragment mass was calculated using 

the experimaitally measured prompt neutron emission data of 
53 

Apalin ^ From the kinetic energies of the fission frag- 

54 their 

ments obtained from Schmitt ^ al. , ^elo cities could be qalcu- 

235 

lated. For thermal neutron induced fission of U , the emission 
of prompt neutrons is almost isotropic. Therefore the emission 
of prompt neutron will not have significant efffect on the velo- 
city of the fission fragment. Hence/ the initial velocity of a 
fission product can be assumed to be the same as that of the 
corresponding precursor fra-ment. These velocities for the 
fission pro ducts ^^Zr and ^^Sr are shown in Table 5.8. 

(b) Effective nuclear charges of fission products - It should 
be noted that a particular fission product of mass n^Jmber A^and 
nuclear charge Z^^can be formed in two ways; 



The initial velocity and the weighted 

average nuclear charge (Z^) for different fission products 
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i) Directly by the emission of prompt neutrons from the 

primairy fission fragment. Pron^t neutron emission from the 

highly excited primary fission fragment occurs in less than 
-13 

*^10 second after the fissile parent nucleus undergoes 
scission. 


ii) By the beta decay of the short-lived fission products 

of mass nxamber A but atomic number The formation 

95 89 

Zr and Sr in this way consists of the following chains: 




8 ^ 140,. 8“^ 140., 6 ^ 140„ 6' 

> Xe — ^ — > Cs ^ Ba — i 


95 

40 


Zr : 


54 


8 95 

» Kr 


36 


55 

8~ 95 , 

^ Rb 

37 


56 

95 


B 95 8 

> '-9r » Y ► 


38 


39 


95 

40 


Zr 


89^ 


8 “ 89 - 

— > Br 

35 


K 89 -, 
— » Kr 




36 


37 


8 89„ 

— 9^ Sr — -> 

38 


Therefore, for any valid comparisons between the calculated and 
the experimentally measured ranges of a fission product of mass 
number A- and atomic number , it is necessary to replace 
in Eq. (5.1) by a weighted average Z^ taking into accomt the 
contribution to its yield from all short-lived fission product 
which undergo beta decay into it. If the independent yield of 
a fission product of mass number A^ and atomic number Z^ is 
designated by Y (Z.), the fractional independeit chain yield 
Ia ( Z^) is given by 
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' ^1 ) 

1-r (Z.)- = - z;- ;/ 

^ ! ^A, (z,.) 


(S.4) 


where the summation includes all fission products having mass 
number A^. The weighted average nuclear charge for the 
fission product of mass number Aj^ is aiven by 


^1 = 


S Z. Ij, (z.) 
Z^ ^ ^1 ^ 


2 I, (Z^) 

Z. 1 
1 


(5.5) 


The fractional independent chain yields were obtained frcan 

55 

Mukherji's prescription. The effective nuclear charges calcu- 
lated in this way for ^^Zr and ^^Zr are shown in Table 5.8* 

140 95 89 

The ranges of the fission products Ba, Zr and Sr 
calculated by using the effective nuclear charges and the 
initial velocities in Eq. (5.1) are given in Table 5.9 and 
are compared with the corresponding experimental values. Since 
Eq, (5.1) for calculating ranges of fission products has be&i 

A 

found to be quite reliable* the good agreement between the 
experimental and calculated values can be taken to be an indirect 
evidence of the relaiability of the present experimental technique, 



Comparison between the experimental 'the 

calculated (R^^i) ranges ^ 
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CHAPTER 6 


SUMMARY 


The ranges of heavy ions in various elemental and conplex 
media have been calculated. For doing this a set of stopping- 
power equations have been developed based on the earlier work 
by Bohr and that by Srivastava and Mukherji, An empirical 
method for calculating the mean ionization potential for all 
elemental media has been given. The "effective quantum nuntoer," 
f(Z) = 0.3634 which has been empirically introduced, 

seems to describe the velocity distributions of electrons in 
atoms vTith good accuracy, as is evident, although indirectly, 
from the good agreement between the calculated and experimental 
stopping-powers and ranges of heavy-ions both for light and 
heavy stopping media. 

140 95 89 

The ranges of the fission products Ba, Zr and Sr 

n O C 

from U-fission have been measured in copper using an electro- 
lytic technique. Though it has been well-known that metals can 
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be cjuantitatively dissolved as an anode of an electrolytic celli> 
it is for the first time that this idea has been used for range 
measurements/ and the method seems to be quite accurate* The 
method can possibly be extended to the case of other metals, 
like Ag, Cr, Pb, Sn and Zr, having low oxidation potentials for 
measuring fission product or heavy ion ranges* 
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APPENDIX 

The chemical corrposition of the "Standard" G-5 emulsion"^ 
has been used for obtaining the empirical formula for G-5 


nuclear 

emulsion of density 

3 

3.815 g/cm . 


Element 

3 

Weight, g /cm , 

of the emulsion 

Number of gram- 

atoms / cm of 
the emulsion 

Number of atoms 

per atoms of Ag 

Ag 

1.8088 

0.01676 

1.000 

Br 

1.3319 

0.01667 

0.995 

I 

0.0119 

0.00095 

0.006 

C 

0.2757 

0.02298 

1.370 

H 

0.0538 

0.05330 

3.185 

N 

0.0737 

0.00525 

0.314 

b 

0.2522 

0.01575 

0.932 

s 

0.0072 

0.00022 

0.013 


The Empirical formula of G-5 nuclear emulsion is/ thus, 
AgBrQ^ggS Iq.OOS *^1.370 ^0.314 ^0.013 '^3.185 ^0.932 * 
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